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Abstract 

The one-loop Minimal Supersymmetric Standard Model (MSSM) contributions 
to the weak and electromagnetic dipole form factors of heavy fermions are reviewed. 
For the Z boson on shell, the weak-magnetic and weak-electric dipole moments of 
the r lepton and the b quark can be defined and directly connected to observables. 
But far from the Z peak, the weak and electromagnetic dipole form factors are not 
enough to account for all the new physics effects. In the context of the calculation 
of the process e + e _ —> tt to one loop in the MSSM, we compare the impact on 
the phenomenology of the CP-violating dipole form factors of the top quark with 
the contribution from CP-violating box graphs. Some exemplificative observables 
are analyzed and the relevance of both the contributions is pointed out. The set 
of tensor integrals employed, the one-loop expressions for the electromagnetic and 
weak dipole form factors in a general renormalizable theory and the SM and MSSM 
couplings and conventions are also given. 


1 Introduction 


The investigation of the electric and magnetic dipole moments of fermions 
provides deep insight in particle theory. The measurement of the intrinsic 
magnetic dipole moment (MDM) of the electron proved the correctness of the 
hypothesis of half-integer spin particles [1] and is one of the most spectacu¬ 
lar achievements of quantum field theory predictions. More precise studies of 
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electron and muon showed afterwards the presence of an anomalous contri¬ 
bution to the MDM (AMDM) and imply very accurate tests of the quantum 
structure of the Standard Model (SM). The measurements of the (g e — 2) and 
( g M — 2) available [2] are in perfect agreement with the SM predictions to sev¬ 
eral orders in the perturbative expansion of the theory (cf. [3] and references 
therein). Furthermore, with the expected precision at the E821 Brookhaven 
experiment [4] it will be possible to improve the previous measurement of 
(< 7 /t — 2) by a factor 20. Therefore the MDMs can be used, together with the 
precision tests at the Z resonance from LEP and SLC and the new results of 
LEP2 and TEVATRON, to set bounds on possible New Physics effects beyond 
the SM [5], 

The importance of the analysis of the electric dipole moment (EDM) of elemen¬ 
tary and composite particles is intimately related to the CP violating character 
of the theory. In the electroweak SM there is only one possible source of CP 
violation, the <5ckm phase of the Cabibbo-Kobayashi-Maskawa (CKM) mixing 
matrix for quarks [6]. Currently the only place where CP violation has been 
measured, the neutral K system, fixes the value of this phase but does not con¬ 
stitute itself a test for the origin of CP violation [7]. On the other hand, if the 
baryon asymmetry of the universe has been dynamically generated, CP must 
be violated. The SM cannot account for the size of the observed asymmetry 
[8]. In extended models (beyond the SM) many other possible explanations of 
CP violation can be given. In particular, in supersymmetric (SUSY) models 
[9,10] CP violation can appear assuming complex soft-SUSY-breaking terms. 
Two physical phases remain in the GUT constrained MSSM [10-12], enough 
to provide the correct size of baryon asymmetry in some range of parameters 
[13]. But the most significant effect of the CP violating phases in the phe¬ 
nomenology is their contribution to the EDMs [14], Unlike the SM, where the 
contribution to the EDM of fermions arises beyond two loops [15], the MSSM 
can give a contribution already at the one-loop level [9]. 

The measurements of the neutron, electron and muon EDMs [16,17] constrain 
the phases and the supersymmetric spectrum in a way that may demand fine 
tuning (supersymmetric CP problem): SUSY particles very heavy (several 
TeV [18]) or phases of 0( 1CT 2 ) [9]. Very large soft-SUSY-breaking masses are 
unappealing as it seems natural to demand the SUSY spectrum to be at the 
electroweak scale.Q On the other side, general universal soft-SUSY-breaking 
terms can be taken only if the genuine SUSY CP phases are vanishing. In this 
case the CP violation is originated via the usual SM CKM mechanism and 
the supersymmetric spectrum affects the observables only through renormal¬ 
ization group equations and/or one-loop radiative contributions. Of course 
in this scenario one has to construct models in which the SUSY phases nat- 


1 Moreover if the SUSY spectrum is in the TeV region this could also give rise to 
relic densities unacceptably large. 
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urally vanish [19] and at the same time provide some other non-standard 
mechanism for explaining electroweak baryogenesis. There exist also ways of 
naturally obtaining small non-zero soft-phases which leave sufficient CP vio¬ 
lation for baryogenesis [12]. But, in general, to get observable effects in most 
electroweak processes one has to relax the assumption of soft-term univer¬ 
sality. Several attempts have been made, following this direction, to use CP 
violation from top-squark mixing: a complex parameter A t would yield large 
CP violating effects in collider processes involving top quarks [20] .Q Besides, 
due to renormalization-group-induced effects on the other low energy phases, 
the phase of A t is constrained by the EDM of the neutron [22], One can also 
make the hypothesis that, due to cancellation among the different components 
of the neutron EDM (constituent quarks and gluons), the SUSY phases can 
still be kept of 0(1) and the SUSY spectrum at the electroweak scale sat¬ 
isfying the experimental bounds [23]. Finally, in [24] it is shown that, in the 
constrained MSSM, large CP violating phases are compatible with the bounds 
on the electron and neutron EDMs as well as with the cosmological relic den¬ 
sities. In view of all these arguments we keep our analysis completely general 
and consider the SUSY CP-phases as free parameters. 

Some attention has also been payed to the study of possible CP violating ef¬ 
fects in the context of R-parity violating models [25]. In this class of models 
new interactions appear providing extra sources of CP violation (still pre¬ 
venting fast proton decay). They can explain the CP violation in the K sys¬ 
tem (with no need of the CKM phase) without introducing anomalous Flavor 
Changing Neutral Current (FCNC) contributions [26]. In any case we assume 
in the following R-parity conservation for simplicity. 

Recently a number of works have been devoted to the analysis of weak dipole 
moments (WDM). The WDMs are defined, in analogy to the usual DMs, 
taking the corresponding on-shell chirality-flipping form factors of the Zff 
effective vertex. The SM one-loop contribution to the anomalous weak mag¬ 
netic dipole moment (AWMDM) has been calculated for the r lepton and the 
b quark in [27,28]. The CP-violating weak electric dipole moment (WEDM) 
is in the SM a tiny three-loop effect. The WDMs are gauge invariant and 
can directly connected to physical observables. While for the r case, using ap¬ 
propriate observables [27,29,30], an experimental analysis is feasible, for the 
b case the situation is complicated by hadronization effects [31]. The SM pre¬ 
dictions are far below the sensitivity reachable at LEP [27] but non-standard 
interactions can enhance these expectations (2HDM [32], MSSM [33]) espe¬ 
cially for the (CP-violating) WEDMs (2HDM, leptoquark models [34], MSSM 
[35]). The experimental detection of non-zero AWMDM or WEDM of heavy 
fermions, at the current sensitivity, would be a clear evidence of new physics 


2 Large non-SM CP violating top-quark couplings could be probed at high energies 
colliders like the NLC [21]. 
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beyond the SM. 


In the perspective of the next generation of linear e + e” colliders, we extend 
the previous analyses on the WDMs to consider the t quark dipole form fac¬ 
tors (DFF) and, in particular, the CP-violating ones. In [36] an independent 
analysis of the t quark EDFF and WEDFF can be found. Since the t quark is 
very heavy one expects this fermion to be the best candidate to have larger 
DFFs. Beyond the Z peak (s > 4mf) other effects are expected to give contri¬ 
butions to the physical observables. In fact the DFFs in a general model are 
not guaranteed to be gauge independent. An exception is the one-loop MSSM 
contribution to the CP-violating electromagnetic and 'weak DFFs (EDFF and 
WEDFF) which do not involve any gauge boson in internal lines. Anyway, 
although this contribution is indeed gauge invariant any CP-odd observable 
will be sensitive to not only the DFFs but to the complete set of one-loop 
CP-violating diagrams involved. In this work we present the full calculation 
of the expectation value of a set of these observables in the context of the 
MSSM and compare the size of the different contributions. 

The paper is organized as follows. In Section 2 we present the general effective 
vertex describing the interaction of on-shell fermions with a neutral vector 
boson. The definitions and the generic expressions of the DFFs for all the 
contributing topologies at the one-loop level are also given. In Section 3 we 
evaluate these expressions for on-shell Z and the case of the r lepton and 
the b quark (their AWMDM and WEDM) in the context of the SM and the 
MSSM. A discussion of the supersymmetric limit is also given. The evaluation 
of the t quark EDFFs and WEDFFs at yfs = 500 GeV is presented in Section 
4. In Section 5 we review observables sensitive to the r AWMDM at the Z 
peak as well as other CP-odd observables valid in general. In Section 6 we 
evaluate specific CP-odd observables for the t quark pair production in e + e~ 
colliders to one loop in the MSSM and compare the influence of the t EDFF 
and WEDFF with that of the CP-violating box diagrams. Our conclusions are 
presented in Section 7. In the Appendices one can End the definition of the 
one-loop tensor integrals and the SM and MSSM couplings and conventions 
that have been used. 


2 The dipole form factors 

2.1 The Vff effective veidex 


The most general effective Lagrangian describing the interaction of a neutral 
vector boson V with two fermions can be written, using at most dimension 
five operators, as a function of ten independent terms: 
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Table 1 

C, P, T properties of the operators in the effective Lagrangian of Eq. (1). Their 
chirality flipping behavior is also displayed. 


Operator 

Coefficient 

P 

CP 

T 

Chirality Flip 


9v 

+ 

+ 

+ 

NO 

y^i|r 7At75 ^ 

9A 

- 

+ 

+ 

NO 

P^-Fi F 

9m 

+ 

+ 

+ 

YES 

P^F 75F 

9e 

- 

- 

- 

YES 

P^Fi d u a IM pV 

9ts 

+ 

- 

- 

YES 

V ^F & J 75 F 

9tp 

- 

+ 

+ 

YES 

(a • P) FF 

9S 

+ 

- 

- 

YES 

(id • P) T75T 

gp 

- 

+ 

+ 

YES 

(^P") ^F 

9TM 

+ 

+ 

+ 

YES 

(iW) *<V75 F 

9te 

- 

- 

- 

YES 


C vff = V»{x)*{x) 


In (flV - 9a 75) + i dp (g M + ifi'ETs) 


+i d u < 7 ^ (g TS + E/tpTs) 


^(x) 


{i& i V v (x))V{x) 


9nv (i<7s + 9p1$) + (ifl'TM + S'teTs) 


*(*)■ ( 1 ) 


The first two coefficients, i.e. gy and g A , are the usual vector and axial-vector 
couplings. They are connected to chirality conserving dimension four opera¬ 
tors. All the other coefficients in Eq. (1) multiply chirality flipping dimension 
five operators and can receive a contribution only through radiative correc¬ 
tions in a renormalizable theory. The operators associated to r/ v , g A , g M , g P , 
c/xm and g^p are even under a CP transformation. The presence of non van¬ 
ishing c/e, PS) fi'TE and gxs yields a contribution to CP-violating observables. 
In Table 1 we summarize the C, P, T and chirality properties of each operator 
introduced in the effective Lagrangian. 

By Fourier transform of Eq. (1) one obtains the most general Lorentz structure 
for the vertex Vff in the momentum space: 


F vff = i 
a 


7m (/v - A 75 ) + {q ~ q)n (/m + I/eTs) + Vn (i/s + Ms) 
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+ (q — <lY a ^ (/tS + i/tp75) + P^CT/iv (i/TM + /te75) ) ( 2 ) 


where q and q are the outgoing momenta of the fermions and p = (q + q) is 
the total incoming momentum of the neutral boson V. The form factors /,; 
are functions of the kinematical invariants. Actually they are more general 
than the coefficients g t . In fact any operator of dimension higher than five 
added to the Lagrangian of Eq. (1) is related to a new coefficient g t . But 
every new coefficient can contribute only to the ten independent form factors 
fi introduced in Eq. (2). The parameters g* and the form factors f % can be 
complex in general. Their real parts account for dispersive effects (CPT-even) 
whereas their imaginary parts are related to absorptive contributions. 

It is possible to lower the number of independent form factors in Eq. (2) 
by imposing on-shell conditions on the fermionic and/or bosonic fields. For 
instance, in the case of on-shell fermions, making use of the Gordon identities: 

2 m/ u 7 M u = ju (q — qY v + u i (q + q) v c t^ v n|, 

2 m/ u 7 m 7 5 v = [u (q + qY 75 v + u i (q - q) u a ^ w 75 uj, 

0 = [u (q + qY v + u i (q - 

0 = {u (q - qY 75 v + u i (q + q )^75 n} , (3) 

one can eliminate /tm, /te, /ts and /tp from the effective Lagrangian. The 
number of relevant form factors can be further reduced taking also the boson 
V on its mass shell. In this case the condition = 0 automatically cancels 
all the contributions coming from fs and fp. The same situation occurs for off- 
shell vector boson V when, in the process e + e _ —» V* —> //, the electron mass 
is neglected. Therefore /g and fp will be ignored in the following. With all these 
assumptions the Vff effective vertex for on-shell fermions is conventionally 
written as: 

r l ff {s) = ie| 7M [V}\s) - A](s) 75 ] + a^q + qY i^^ - ^^75 |, (4) 

where e and m/ are respectively the electric unit charge and the mass of the 
external fermion. The form factors in Eq. (4) depend only on s. As men¬ 
tioned above Vj'f s) and A / (s) parameterize the vector and axial-vector cur¬ 
rent. They are connected to the chirality conserving CP-even sector. The 
form factors ct/ (s) and d / (s) are known respectively as anomalous magnetic 
dipole form factor (AMDFF) and electric dipole form factor (EDFF). They 
are both connected to chirality flipping operators. In a renormalizable theory 
they can receive contribution exclusively by quantum corrections. The EDFFs 
contribute to the CP-odd sector and constitute a source of CP-violation. 
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Class I: Class II: Class III: 



Fig. 1. The one-loop Vff diagrams with general couplings. 

The dipole moments (DM) are defined taking the corresponding vector bosons 
on shell, s = My. For V = 7 one gets the usual definitions of the pho¬ 
ton anomalous magnetic dipole moment (AMDM) and electric dipole moment 
(EDM). The definitions of electric charge, magnetic and electric dipole mo¬ 
ments [37] consistent with our convention for the covariant derivative (B.l) 
are respectively: 


charge = —e Vj(0) — e Qf , (5) 

MDM=^(V7(0) + aJ(0)), (6) 

EDM = dj(0) . (7) 

Thus the AMDM of a fermion is aj(0) = — Qf(g / — 2)/2 with gf being the gy- 
romagnetic ratio. The axial-vector coupling AJ vanishes. For V — Z, the quan¬ 
tities aj = af (Mf) and df = dj (Mf) define the anomalous weak-magnetic 
dipole moment (AWMDM) and the weak-electric dipole moment (WEDM ). 3 


2.2 One-loop generic expressions of the dipole form factors 

All the possible one-loop contributions to the Oj (s) and dj (s) form factors 
can be classified in terms of the six classes of triangle diagrams depicted in 
Fig. 1. The vertices involved are labelled by generic couplings for vector bosons 
yd) — W tJ , l-'l-'d, fermions and scalar bosons t&fc, according to the 

following interaction Lagrangian: 


C = ieJ(Wl u W^V u - W^Wjy, + V^WlW u ) + i eGjkV^j d„ 

3 Massless and neutral fermions may have magnetic moments. The usual 
parametrization in (4) must be generalized by the replacement e aj (s)/2mj —> 
fi V (s). The dipole moments of neutrinos are then given by pX (My). 
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+ {e* f (S jk - P jk ^k^j + eK jk V^ k % + h.c.} 

+eVi k %'f(Vjt ) ~ 4^ 6 )tf, . 


(8) 


Every class of diagrams is calculated analytically and expressed in terms of 
the couplings introduced in (8) and the one-loop three-point integrals C t (see 
App. A). The result is given in the’t Hooft-Feynman gauge. 

• [Class I]: vector-boson exchange: 



In the case of gluon exchange one has to substitute a for a s and vj^ for Ti 
(the SU(3) generators) and take = 0. The sum over the index / yields 
a color factor Cf = 4/3. The gluon contribution to the CP-violating form 
factor dj (I) vanishes in general. 

• [Class II]: fermion exchange and two internal vector bosons: 



[Class III]: scalar exchange: 



Q/(s) 

2m f 


m = Tz{2m,Y.^ 

jkl 


4 W(m )= f 

e 47 t 


VjP(SijS; t + p^) 

+A% ) (S,,P,l + PyS,;)] [2C+ - c* 
-2^m t Re[vP(S l: iS; k - P tj P? t ) 

-A%\s,jP: t - PijS; k )] [ct + ct] kj } 
vP\PiA + s,,p; k ) 


jkl 


J kjl 


(13) 


| — 2 m f ^ Im 

jkl 


+41 \s v s; k + p v pt k )] her - c r 


kjl 


+2'£,m k lm[v}P(P lj S; k - S„P? k ) 

jkl 

+A%>(s, j si, - fljfjy] [c, + + cr] w } 

[Class IV]: fermion exchange and two internal scalars: 


(14) 


Q-/ 0 ) 

2 mi 


/Y/ r _ r -i 

(iv) = -ld 2m ' ?ad G i‘( s i lS « + ^P5)] [2C+ - cf 


jkl 


J kjl 


+ '£m l Re[c jk (S jl Sl l - PjiPt,)} [2 Ct - C„] w } (15) 


jkl 


®V)=f 


12 m f ^2 I m 

jkl 


GASflPZ + PflSru)} \2Ct-c-, 


5>,Im - PflSi,)} [2C+ - C„] } (16) 

jkl 


J kjl 
kjli 


[Class V+VI]: fermion exchange, one vector and one scalar internal boson: 
aX(s) 


V 

2 m f 

d V f (s ) 


a 


(V + VI) = £2£Re[A>(V«S*, + A^P’M [C+ + Cfl (17) 


(*); 


47T 


jkl 


J kjl 


(v + VI) = -^ 25 im[jr jt (v<‘>p;, + 4 ?sj)] [c+ + c r ] jis) 


jkl 


kjl 


In Eqs. (9-18) the shorthand notation [C]kji stands for the three-point tensor 
integrals C(—q,q,Mk,Mj,Mi). The integrals appearing in the previous Eqs. 
are UV and IR finite. All the expressions are proportional to some positive 
power of a fermion mass, either internal or external, consistently with the 
chirality flipping character of the dipole moments.^ For class V and VI dia¬ 
grams the mass dependence is hidden in the product of the Yukawa couplings 

1 This is as expected when applying the mass-insertion method: to induce a flip 
in the fermion chirality one introduces a mass in either one of the internal fermion 
lines, picking a mass term from the propagator, or in the external fermion lines, 
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Fig. 2. SM diagrams for the Zbb vertex (Higgs sector excluded). 

Sij ( Pij ) and the dimensionful parameter K t] . Hence, the heaviest fermions 
are the most promising candidates to have larger DFFs. Eqs.(9-18) also show 
that, in general, the DFFs for massless fermions are not vanishing but propor¬ 
tional to masses of fermions running in the loop. The SM cancellation of the 
massless neutrino DFFs is only due to the absence of right-handed neutrinos. 
Finally, notice that all the contributions to the EDFFs are proportional to 
the imaginary part of certain combinations of couplings. A theory with real 
couplings has manifestly vanishing EDFFs. 


3 SM and MSSM predictions for the r and b WDMs 


In the previous Section we derived the DFFs for the generic interaction La- 
grangian of Eq. (8). In this Section we focus in the calculation of the WDMs in 
two specific models: the SM and the MSSM. These quantities, defined at the 
Z peak, are gauge invariant and can be directly related to observables. The 
r lepton and the b quark are the heaviest fermions to which a Z boson can 
decay and hence the optimal candidates to have larger weak dipole moments. 


3.1 The t and b magnetic and electric WDMs in the SM 


In the clectroweak sector of the SM, adopting the ’t Hooft-Feynman gauge, 
there are 14 diagrams contributing at one loop to the Zff boson interaction. 
One more diagram has to be considered (the gluon exchange) if also QCD 
is included. All the possible SM diagrams for the Zbb vertex are given in 
Figs. 2 and 3. The calculation in Ref. [33] is in perfect agreement with the 
results of Refs. [27,28] for both the r lepton and the b quark. Taking as input 
parameters m T = 1.777 GeV, = 4.5 GeV, rn t = 175 GeV, M z = 91.19 GeV, 
s'^y = 0.232, a = 1/128 and a s = 0.118, the pure electroweak contributions 
are: 


a”[EW] = (2.10 + 0.61 i) x 10 -6 , (19) 

a™ [EW] = [(1.1; 2.0; 2.4) - 0.2 i] x 10~ 6 , (20) 

using the equations of motion for the free fermion. 
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for three different values of the Higgs mass (respectively Mh o = Mz, 2Mz 
and 3 Mz)- Note that for the r case the Higgs boson contribution has a very 
small impact, changing the final result by less than 1 % in the range of the 
Higgs boson masses considered. A much stronger dependence is observed for 
the b case. The CKM matrix element V tb is taken equal to one and all the off- 
diagonal entries vanishing. The whole SM contribution (EW+QCD) for the b 
quark is dominated by the gluon exchange: 

a™ [SM] = (-2.98 + 1.56 i) x 1(T 4 . (21) 


This large enhancement is as expected by the replacement a —> a s . 

The only phase present in the SM, the <5ckm, is not sufficient for generating 
one-loop contributions to the (W)EDM.[[] One needs to go beyond the two- 
loop level [15]. A very crude estimate of the tiny three-loop SM contribution 
to the (W)EDM can be done using simple power counting arguments. For a 
general fermion with mass rrif dj[SM] is of the order of eG^m/a 2 ^ J/(47t) 5 
with J = C 1 C 2 C 3 S 152533,5 being an invariant under reparameterizations of the 
CKM matrix [14]. For the r lepton and the b quark we have: 

d™ [SM] < 1.5 x 10 “ 19 ^ and d^SM] < 1.5 x 1(T 18 fi b . (22) 

We use the scale unit /aj = e/2 nif (“magneton”) with /i T (n b ) = 5.55 (2.19) x 
ICC 15 ecm, respectively. 

3.2 The r and b AWMDMs in the MSSM 


The particle content of the MSSM includes the SM spectrum, the SUSY part¬ 
ners and an extended Higgs sector consisting of two doublets. We assume 

5 We implicitly assume here a vanishing $qcd phase. 
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Fig. 4. MSSM Higgs diagrams for Zbb. 



implicitly R-parity conservation. The new non SM diagrams are depicted in 
Figs. 4 and 5. They will be grouped into sets of diagrams including Higgs 
bosons, neutralinos, charginos and gluinos, as in Ref. [33]. All of them belong 
to the topologies of classes III, IV, V and VI. 

The MSSM involves a large set of new parameters currently unknown but 
constrained by present experimental data. For sake of simplicity we assume 
the usual GUT relations for the soft-breaking gaugino mass terms, 

i CM 

M 1 = - tan 2 9 w M 2 , M 3 = —s 2 w M 2 (23) 

3 a 


and a common scalar fermion mass parameter. Moreover we restrict ourselves 
to two typical scenarios, low and high tan/?, set respectively to tan (3 = 1.6 
and tan/? = 50.[^]We also assume generation-diagonal, not universal, trilinear 
soft-breaking terms in order to prevent large FCNC.Q We are left, for given 
values of tan /?, with the following free parameters: the gaugino mass parameter 
for the SU(2) sector M 2 , the Higgs-higgsino mass parameter /i, the mass of the 
pseudoscalar Higgs boson M 4 , the common scalar lepton soft-breaking mass 

6 Currently our adopted value for the low tan /? scenario is close to be excluded by 
Higgs boson searches [41] but we keep it for reference. 

7 For a review about FCNC effects see for instance [26]. 
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parameter mp the common scalar quark soft-breaking mass parameter rriq 
and the trilinear soft-breaking terms Af (for / = b, t or r). The conventions 
for couplings and mixing are extensively explained in the Appendices. In this 
Section all the supersymmetric parameters are taken real. 


3.2.1 MSSM Higgs contribution to the AWMDMs 

The MSSM Zbb diagrams involving Higgs bosons are shown in Fig. 4. There 
G° and G ± are the would-be-Goldstone bosons and h, H, H ± and A are the 
physical scalar and pseudoscalar MSSM Higgs bosons. Actually not all these 
diagrams contribute to the AWMDM: the class IV diagrams with neutral Higgs 
bosons (necessarily of opposite parity) identically vanish, also in the SM, due 
to the way their couplings to fermions (B.78-B.83,B.85) combine in (15). There 
are only three independent parameters involved in the Higgs sector: tan /3, the 
mass of the pseudoscalar Ma, and the common scalar fermion mass parameter 
mj. The Higgs contribution decreases with growing Ma and mj, consistently 
with the decoupling theorem [38,39]. 

The Higgs contribution to the Re (a"’) coming from the diagrams of classes 
III and IV are proportional to (m T /Mz) 4 tan 2 (3 and those from the diagrams 
of classes V and VI are proportional to (m T /Mz) 2 . Due to the fourth power 
suppression the sum is not very sensitive to tan (3. We obtain 

Re (a™ [Higgs]) = —0.3 (—0.4) x 10~ 6 for tan/3 = 1.6 (50), (24) 


mq = 250 GeV and Ma = 100 GeV. Assuming the present lower bounds for 
the masses of the Higgs bosons [40,41] the only contribution to the Im(«"') 
can arise from class III diagrams and it is in general smaller than the SM one. 

The same considerations apply for the Higgs contribution to the b AWMDM. 
The only relevant difference is the presence of a dominant term proportional 
to ( mbirit/Mz ) 2 . The Re(a™) is: 

Re (a™[Higgs]) = —3.8 (+0.8) x 10 -6 for tan/3 = 1.6 (50), (25) 


rriq = 250 GeV and Ma = 100 GeV. Again the contribution to the Im(a^) is 
negligible. 


3.2.2 Chargino contribution to the AWMDMs 

There are two diagrams with charginos and scalar fermions (Fig. 5). In the re¬ 
gion of MSSM parameters not ruled out by present experiments [41] charginos 
and scalar fermions cannot be pair produced in Z decays. Therefore they can 
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(b) 


Fig. 6. Re (a™ [char]) for low (a) and high (b) tan [3 in units of 10~ 6 in the plane 
M 2 — p. The contour solid-lines correspond to the lightest chargino masses m^± = 90 
and 200 GeV and the contour dashed-lines correspond to the lightest neutralino 
masses m^o = 15 and 100 GeV. The common scalar quark mass parameter is fixed 
to mq = 250 GeV and m t LR = 0. From [33]. 


contribute only to the real part of the AWMDM. The free parameters involved 
are tan /3, M2 , / 1 , mj (and A t in the b case, since there is no mixing in the scalar 
neutrino sector). Increasing the scalar fermion common mass only produces 
decoupling of the scalar fermion from the spectrum. In the large M 2 and \n\ 
region also the charginos decouple. 

The chargino contribution to the real part of the r AWMDM is: 

Re (df [char]) = ±0.2 (±7.0) x 10 ' 6 for tan/3 = 1.6 (50), (26) 


with rriq = 250 GeV and M 2 = \/a\ = 200 GeV. The enhancement for higher 
tan (3 is due to the dependence of the chargino and neutralino couplings on 
the Higgs vacuum expectation values. 

For the real part of the b AWMDM one gets: 

Re(a™[char]) = ±1.1 (±33) x 10 ~ 6 for tan/3 = 1.6 (50), (27) 


with rriq = 250 GeV, M 2 — \pi\ — 200 GeV and A t = /icot/3. The ± stands 
for sign(/i). Varying the value of A t in the range ±/i cot /3 there is only a small 
effect in the AWMDM. We thus choose the value of A t that makes the off- 
diagonal entry of the t scalar quarks mixing mass matrix vanish. In such a 
way the scalar fermion mass eigenstates are always physical. Fig. 6 shows the 
contour plot of the Re(a™[char]) in the M 2 —/1 plane. The contour lines for fixed 
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masses of the lightest chargino (90 and 200 GeV) and lightest neutralino (15 
and 100 GeV) are also given for illustration. The contributions are enhanced 
by increasing tan /3. In the low tan /3 scenario the contribution is of order 10” 6 
(Fig. 6a), while in the high tan (3 scenario it becomes of order 10~ 5 (Fig. 6b). 
The contour plots for the r AWMDM are similar to the b quark ones but 
suppressed by factors between m T /mb and (m T /rrib) 2 [see Eqs. (13) and (15)]. 


3.2.3 Neutralino contribution to the AWMDMs 

The two diagrams involving neutralinos and scalar fermions are shown in 
Fig. 5. The scalar fermion masses are bounded by present experiments to 
be heavy but the neutralino masses can be lighter than half the mass of the 
Z [41]. This allows for the possibility of a contribution to the imaginary part 
of the AWMDM through class III diagrams. A non negligible effect can arise 
only extremely near to the neutralino production threshold. The neutralino 
contribution to the AWMDM depends on tan/3, M 2 , /i, rrij, as well as the 
trilinear soft-breaking mass term Af. No sizeable deviations due to changes 
in Af in the range ±/i tan /3 are observed. To avoid an unphysical mass for 
the fermion we take m RR = 0. The usual decoupling properties for scalar and 
fermions have been successfully checked. 

For the real parts of the r and b AWMDM one gets: 


Re «[neut]) = -0.02 (-0.4) x 1(T 6 , (28) 

Re (a™[neut]) = — 0.2 (—10.5) x 10“ 6 for tan /3 — 1.6 (50). (29) 

and M 2 = ft = 200 GeV. Like in the chargino contributions the result is 
enhanced with increasing tan /3. 


3.2.4 Gluino contribution to the b AWMDM 

In the case of the b AWMDM one more diagram must be considered: the one 
with a gluino and two b scalar quarks running in the loop. Due to the a — » a s 
replacement, we expect from the gluino sector a large contribution to the b 
AWMDM. The present bounds on b scalar quark masses [41] only allow a real 
contribution to the AWMDM. The free parameters are tan/3, niq, m,g , A& and 
/i. The last two affect the result only through the off-diagonal term of the b 
scalar quark mass matrix m b LR = Ab — /i tan /3. This parameter is responsible 
for the mixing in the b sector and intervenes in association with the chirality 
flipping in the gluino internal line. Thus the contribution to the AWMDM is 
almost proportional to m b LR . We assume /i tan fj as typical scale for m b LR . 

In Fig. 7 the gluino contribution to the b AWMDM is displayed in the m b LR —mg 
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Fig. 7. Re(af [glu]) low (a) and high (h) tan/3 in units of 10~ 6 in the plane m b LR —mg. 
The common scalar quark mass parameter is fixed to nig = 250 GeV. The region 
above the upper horizontal line is unphysical (m~ <0). From [33]. 


plane. In the low tan/3 region (Fig. 7a) the value of m b LR does not affect 
significantly the mass of the lightest b scalar quark (always above 200 GeV). 
For zero gluino mass, only the term proportional to the mass of the b quark 
provides a contribution. As we increase the gluino mass, the term proportional 
to rrig dominates, especially for large m b LR , becoming again suppressed at high 
nig due to the gluino decoupling. For high tan /3 the behavior is analogous but 
larger values for the AWMDM can be obtained as shown in Fig. 7b (extremely 
high values of m b LR are excluded due to unphysical m~ < 0). Typical values 
are: 


Re (a™[glu]) = —1.1 x 10 6 for tan/3 = 1.6, m b LR = 200 GeV, (30) 
Re (a™ [glu]) = —18.9 x 10” 6 for tan/3 = 50, m b LR = 7 TeV, (31) 

with the gluino mass fixed to rrig = 200 GeV. Using the GUT constraint (23) 
higher values for the gluino mass ( rrig = M 3 ) seem more likely and in this case 
the results in Eqs. (30, 31) are approximately reduced linearly with trig. 


3.2.5 MSSM total contribution to the AWMDMs 

The imaginary parts of the r lepton and b quark AWMDM are provided only 
by the diagrams with Higgs bosons and neutralinos. The total contribution in 
a reasonable region of the MSSM parameter can reach at maximum 10” 6 and 
is typically lower than the corresponding SM contribution. 

The real part is dominated by the chargino and the gluino sectors, when 
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present. It is roughly proportional to tan/3. The total contribution can reach 
the values: 


|Re«[MSSM])|~0.5 (7) x 1(T 6 , 
|Re(a™[MSSM])| ~2 (50) x RT 6 


(32) 

(33) 


for tan/3 = 1.6 (50) respectively. In the low tan/3 scenario the MSSM contri¬ 
bution to the AWMDM is of the same order than in the SM (19,20). More 
interesting is the high tan/3 scenario where the MSSM contribution to the 
real part can reach values one order of magnitude higher than the pure elec- 
troweak SM contribution, especially in the b quark case (20), but still a factor 
five below the standard QCD contribution (21). 


3.3 The t and b WEDMs in the MSSM 

Eqs. (9-18) show that for having non-vanishing contribution to the WEDMs 
one has to deal with complex couplings. Many of the parameters in the SUSY 
Lagrangian can be complex but not all the phases originated are physical. 
Some of them can be absorbed by a redefinition of the MSSM fields (see App. B 
for a complete discussion of the procedure). Our choice of CP-violating phases 
is:Q (pp = arg(//), ipj = arg(m I / R ) (/ = r, t, b ) with m^ R = A t - fi* cot/3 and 

uq/ R = A Tib - /A tan/3. 

Since the MSSM contains a CP-conserving Higgs sector'/] and the SM provides 
contributions to the (W)EDMs beyond two loops, the relevant diagrams are 
the genuine SUSY graphs of classes III and IV [35]. 

3.3.1 Chargino contribution to the WEDM 

There is only one phase, involved in the chargino contribution to the r 
WEDM, as there is no mixing in the scalar neutrino sector. The result grows 
with tan f3. It also depends on the common scalar lepton mass (whose effect 
consists of diminishing the result through the tensor integrals) and the /i and 
M 2 mass parameters. Taking M 2 = \/a\ = m z ~ = 250 GeV one obtains: 


Re (d"’[char]) = 0.18 (5.52) x 10 6 fi T for tan/3 = 1.6 (50) (34) 


8 Such a choice leads to a dependence on of chargino and neutralino masses. 
Conversely the scalar fermion masses are independent of <pj. 

9 A one-loop non vanishing contribution to the WEDM is possible in a general 
2HDM [42], 


17 




Fig. 8. Re(df [char]) for low tan (5 (a) and Re(d™[neut]) for high tan (3 (b) in the 
plane <pj — in units of 1CT 6 pj, with M 2 = \p\ = rriq = 250 GeV. The mixing 
\ mt Ln\ = |//1 cot/? and \m b LR \ = |/i|tan/3, respectively. From [35]. 


in the limit = tt/2, for which the WEDM takes the largest value. 

Two CP-violating phases are involved in the b case: tp^ and p R In Fig. 8 a 
the dependence on these phases is shown, for M 2 — \p\ — rriq = 250 GeV, 
= |/i cot /? and both low and high tan /? scenarios. The maximum effect 
on the WEDM is obtained for = tt/2 and pi = it. One gets up to 

Re (eschar]) = 1.17 (27.1) x 10 -6 /xj, for tan/3 = 1.6 (50). (35) 


In the high tan/3 scenario our assumed takes a small value and the 

dependence on iff tends to disappear. To have an idea of the maximum value 
achievable for the chargino contribution, we show in Fig. 9 the dependence 
on M 2 and /i| for <p /t = ±7t/2 and m t LR = 0 with the previous value for the 
common scalar quark mass parameter. Assuming the present bounds on the 
chargino masses [41], there is no contribution to the imaginary part of the r 
and b WEDM. 


3.3.2 Neutralino contribution to the WED Ms 

Now both and (p f contribute. Assuming that \m T LR \ is of the order of 
\fi | tan/3 or below, there is no large influence of ip? on the neutralino contri¬ 
bution, for both low and high tan /3. This contribution is roughly proportional 
to sin (fin regardless of the value of (pf. Taking ip^ = tt/2, M 2 = |//| = 250 GeV 
and mi = 250 GeV one finds: 

Re (d™ [neut]) = —0.01 (—0.25) x 10 ~ 6 /i T for tan/3 = 1.6 (50). (36) 
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Fig. 9. Re (d b [char]) for low (a) and high (b) tan (3 in the plane M 2 —Im(n), in 
units of 10 -6 /if, with |siny> M | = 1 and rn,q = 250 GeV. The mixing |to£ R | = 0. 
The contour solid-lines correspond to the lightest chargino masses m-± = 80 and 
200 GeV and the contour dashed-lines correspond to the lightest neutralino masses 
rri^o = 15 and 100 GeV. From [35]. 


The two relevant CP-violating phases for the neutralino contribution to the b 
WEDM are: p M and p~ h . As before, the most important effect from p~ h arises 
when the off-diagonal term is larger, which in this case corresponds to high 
tan /?, as the trilinear soft breaking parameter is taken to be of the order 
of |/i|ta a/3. The maximum value for the neutralino contribution occurs for 
Tu = Tb = 77 (Fig- 8 b). The total contribution increases with tan/?. Thus 
taking M 2 = |/i| = rriq = 250 GeV and \m b LR \ — |/i tan/? one gets: 

Re (d™[neut]) = —0.29 (—12.6) x 1CT 6 /p, for tan/? = 1.6 (50). (37) 


For presently non excluded masses of the neutralinos [41] there can be a con¬ 
tribution to the imaginary part of the order of 10” 6 (if. 


3.3.3 Gluino contribution to the b WEDM 

The gluino contribution is affected only by m b LR , rriq and the gaugino mass 
nig. Therefore the maximum value occurs for p~ b = 7t/2. The mixing in the b 
sector is determined by m b LR and intervenes in the contribution due to chiral¬ 
ity flipping in the gluino internal line (the contribution proportional to mg). 
The contribution to the WEDM is enhanced by the largest values of \m b LR \ 
compatible with an experimentally not excluded mass for the lightest b scalar 
quark. For zero gluino mass, only the term proportional to the mass of the b 
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quark provides a contribution. Thus for <pi = vr/2 one gets: 


Re (d™[gluinos]) = 0.26 (9.31) x 10 6 pb for tan/3 = 1.6 (50), (38) 


\m b LR \ = M tan/3, M 2 = | p | = rriq = 250 GeV and rrig = M 3 fulfilling the 
GUT relation (23). There exists the possibility for a gluino contribution to 
the imaginary part of the b WEDM if the hypothesis of light gluino is not 
discarded. 


3.3.4 MSSM total contribution to the WEDMs 

The Higgs sector does not contribute and chargino diagrams are more im¬ 
portant than neutralino ones. Gluinos are also involved in the b case and 
compete in importance with charginos. In the most favorable configuration 
of CP-violating phases and for values of the rest of the parameters still not 
excluded by experiments, these WEDMs can be as much as twelve orders of 
magnitude larger than the SM estimates (22): 

|Re«)| < 0.2 (6) x 10” 6 p T (39) 

|Re(d™)| < 2 (35) x 10 -6 pb (40) 

There may be a contribution to the imaginary part if the neutralinos are light 
but this contribution is at least one order of magnitude smaller than the real 
part of the r or b WEDM. 


3-4 Cancellation of the dipole moments in the supersymmetric limit 


A general Vff interaction is restricted to the form (2) by Lorentz invariance. 
Since the Lorentz algebra is a subalgebra of the supersymmetry algebra, this 
interaction is even more constrained in a theory with unbroken supersymme¬ 
try. In fact, supersymmetric sum rules are derived that relate the electric and 
magnetic multipole moments of any irreducible N — 1 supermultiplet [43]. 
Applied to the Vff interaction between a vector boson coupling to a con¬ 
served current and the fermionic component of a chiral multiplet, these sum 
rules force the gyromagnetic ratio to be gf = 2 and forbid an electric dipole 
moment: 


Of = dy — 0 . 


(41) 


The Lagrangian of the MSSM is supersymmetric when the soft-breaking terms 
are removed. For non zero value of the p parameter the Higgs potential has 
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only a trivial minimum. Therefore to keep the particles massive and super- 
symmetry preserved at the same time the choice n — 0 is necessary. Then 
the Higgs potential is positive semi-definite and it has degenerate minima 
corresponding to v = V\ = v 2 (tan/? = 1). The value of Ma = 0 follows 
from such a configuration. Finally the value of the common v is fixed by the 
phenomenology: the muon decay constant and the gauge boson masses. 

In this supersymmetric limit the above mentioned sum rules are valid and 
the magnetic and electric dipole form factors have to cancel. To verify our 
expressions we checked this for the AWMDM of the b quark [35]. Choosing 
the parameters Ab = A t = M 2 = M 3 = 0, /j, — 0, tan (3 = 1 and Ma = 0, 
the SM gauge boson contribution to af [28] is indeed cancelled by the MSSM 
correction including the two Higgs doublets: the gluon and gluino contribution 
cancel among themselves and the neutralinos and charginos cancel the gauge 
boson and Higgs contributions. 


4 SM and MSSM predictions for the top quark dipole form factors 

4.1 SM 


The electroweak contributions to the magnetic and weak-magnetic dipole form 
factors for off-shell gauge bosons are gauge dependent. The pinch technique 
[44] could be used to construct gauge-parameter independent magnetic dipoles 
in the class of R g gauges [45,32] but the prescription is not unique [46] and 
these quantities cannot be observable by themselves. The QCD contributions 
(gluon exchange) to the t (W)MDFF are gauge independent and comparable in 
size to the electroweak predictions at yfs = 500 GeV in the’t Hooft-Feynman 
gauge due to the large mass of the t quark [32], 

There is no contribution to the electric and weak-electric dipole form factors 
to one loop in the SM. 


4.2 MSSM 


The triangle diagrams with SUSY particles are gauge independent by them¬ 
selves (no gauge or Goldstone bosons involved) but the ones including Higgs 
scalars and gauge bosons in the loop are not sufficient to keep the gauge invari¬ 
ance in the case of the magnetic dipole form factor. The CP-violating dipole 
form factors (electric and weak-electric) for which the Higgs sector of the 
MSSM is irrelevant, on the other hand, can be regarded as gauge independent 
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quantities. These form factors have also been considered in [36] 

We investigate the SUSY contributions to the t electric and weak-electric 
dipole form factors performing a parameter scan for which a fixed value y/s = 
500 GeV has been chosen. □ It is important to point out that the region of 
the supersymmetric parameter space that provides the maximal contributions 
varies with s due to threshold effects. We scan the mass parameters M 2 and /i 
in a broad range and the CP-violating phases: pi and p~ h (see Figs. 10 and 

11). The gluino mass is given by the GUT constraint (23). We adopt a fixed 
value for the common scalar quark mass rriq = 200 GeV: this is a plausible 
intermediate value; larger values decrease the effects. Finally, the moduli of 
the off-diagonal terms in the t and b mass matrices are also chosen at fixed 
values \m? LR \ = \m b LR \ = 200 GeV, to reduce the number of free parameters. 
The results are expressed in t magnetons = e/2 m t = 5.64 x lO^ 1 ' ecm. 

The MSSM contributions to the top quark EDFF and WEDFF include: 


4-2.1 Neutralinos and t scalar quarks 

They provide typically small contributions but quite sensitive to the value of 
both the phases involved, p p , and pi (Fig. 10). 

The results are larger for low tan (3 since the chirality flipping mass terms are 
dominated by the t quark, yielding a term proportional to m t cot (3. A term 
proportional to the neutralino masses is also present as well as a negligible one 
proportional to tan (3. The contributing diagrams belong to the classes III 
and IV for the Z case and only to class IV for the 7 case, as the neutralinos 
do not couple to photons. 

As a reference we take the representative values M 2 = p = 200 GeV and 
p p = —pi = 7 t /2. For these inputs the results are 

d/[x°] = ( 0.080 + 0.081 i) x 10“ 3 p t (42) 

df[x°) = (-0.324 + 0.223 i) x 10" 3 p t (43) 

4-2.2 Charginos and b scalar quarks 

As in the case of the neutralinos, the influence of chargino diagrams is en¬ 
hanced for low tan (3. 

10 The results in [36] after revision are in agreement with the ones presented here 

[ 47 ]- 

11 We use the running coupling constants evaluated at y/s = 500 GeV, a s = 0.092, 
a = 1/126. 
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Fig. 10. Neutralino contribution to the real and imaginary parts of the t EDFF and 
WEFF [in 10~ 3 /it units] in the plane for tan/? = 1.6 and the reference 

values M 2 = \p\ = |m^ R | = 200 GeF at \/i = 500 GeV. 

The results depend very little on <p~ b and mainly on the maxima being close 
to ipn = ±7t/2. For increasing tan/? the dependence on ip~ h grows, as it comes 
with a factor proportional to tan/3. 

The chargino contributions are the most important ones. In Fig. 11 the de¬ 
pendence on Im(/i) and M 2 is displayed for tan/3 = 1.6. The only relevant 
CP-violating phase here has been set to the most favorable case, = n/2 
(the negative values of Im(/i) correspond to = —7 t/2 ). The symmetry with 
respect to Im(//) = 0 in Fig. 11 reflects the approximate independence of p b , 
here set to n/2. □ The same does not happen for the neutralino contributions 

12 All the contributions flip sign when the set {p^^p^pf) is rotated by n. They 
vanish accordingly when all the phases are zero. See Fig. 10 for illustration. 
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Fig. 11. Chargino contribution to the real and imaginary parts of the t EDFF and 
WEFF [in 10 _3 /p units] in the plane M2 — Im(/x) for tan/l = 1.6, \m b LR \ = 200 
GeV and ip^ = tt/2 at s/s = 500 GeV. The lower (upper) solid isolines correspond 
to m-± = 90 (200) GeV and the dashed isolines to rn.~o = 25 (100) GeV. 


for a fixed value of (pi ^ 0,7r in the plane M 2 — Im(/i). The plots of Fig. 11 
exhibit a tendency to decoupling of the supersymmetric effects for increasing 
values of the mass parameters. The isolincs for a couple of masses of the light¬ 
est charginos and neutralinos in the same plane are also given for orientation. 
The current LEP2 experimental lower limits are m^o > 25 GeV and m-± > 90 
GeV [41], 

The chargino contributions for the values M 2 = \/i\ = 200 GeV and ip^ = 7t/2 
are 


dj[x ± ] = (0.869 - 1.870 i) x 10~ 3 /p 


(44) 
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ilf[C] = (0.793 - 2.524 i) x 1(T 3 />., 


(45) 


4-2.3 Gluinos and t scalar quarks 

Their effect is roughly proportional to |m^ R |sin iff times a chirality flipping 
fermion mass, either m t or M 3 (the glnino mass). It is damped for heavy gluinos 
circulating in the loop and also for large scalar quark masses (decoupling). 
Both terms have opposite sign to Im(m^ R ) and the one proportional to the 
gluino mass dominates. 

The result for M 2 = 200 GeV and yy = — 7 t /2 is 

(T t [g\ = (0.457 + 0.170 i) x 10 ~ 3 /.q 
df[g\ = (0.155 + 0.059 i) x 10 ~ 3 /.q 

4-2-4 Total contribution 

In view of these results, we establish a set of SUSY parameters and phases for 
which nearly all the contributions sum up constructively at a/s = 500 GeV. 
Our choice is 


(46) 

(47) 


Reference Set #1 : tan/3 = 1.6 

M 2 = \n\ = rriq = |m// = \m b LR \ = 200 GeV 

Vr = -Vi = ~Vl = tt/2, (48) 

for which the t EDFF and WEFF reach the values: 

d? t = (1.407 — 1.618 i) x 10 -3 /j, t (49) 

df = (0.624 - 2.242 i) x 10" 3 q f , (50) 

at a/s = 500 GeV. To illustrate how the maximum effects appear we display 
in Fig. 12 the individual contributions as a function of /s. The masses of 
the supersymmetric partners in the loops are such that there are threshold 
enhancements in the vicinity of /s = 500 GeV. 


5 Observables at the Z resonance: the r lepton case 


We analyze below a set of observables sensitive to the weak dipole moments at 
the Z peak. In this environment the effects of electromagnetic dipole moments 
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Fig. 12. The different contributions to the t (W)EDFF [in 10 3 //t units] for the 
reference set of SUSY parameters of Eq. (48). 


to the process e + e —> // are negligible, as the pair production proceeds 
through the Z resonance. The 7 —exchange and box diagrams are suppressed. 


5.1 The Z width 


The Z width is given in terms of generic vertex form factors by 


T(Z->//) = 


G F Ml 


A 


1 - 


4 m 2 f 

Jq 


1 + 


2 mV 


6 ^ 

-6vfSwCwRe(af) + 4s^/C^/|aJ| 2 


+ Clf I 1 — 


Ami 

Mf 


'Ml 

^8 rrij 


+ 1 
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( 51 ) 


+2s 2 2 I fjw 12 / M Z _ A 

where we have introduced the dimensionless parameter df = 2 mfdj/e for 
simplicity. 

From the previous expression one observes that the real part of the AWMDM 
contributes linearly to the total width (although it is not necessarily the dom¬ 
inant term in the r case for af near the experimental upper limit, due to 
the small value of v r and the M^/mf enhancement factor). The WEDM con¬ 
tributes only quadratically to the total Z width. 

One might get some upper bounds on Re (a™) and \df\ from the Z partial 
widths measurements at LEP and SLC using the previous expressions. Actu¬ 
ally this is not the most appropriate way to put limits to the dipoles as much 
information is integrated out. Moreover the weak-magnetic and weak-electric 
contributions cannot be disentangled and some other effects can interfere the 
measurements. In contrast, the differential cross section for the production of 
polarized fermion pairs contains all the ingredients needed and is the starting 
point for the construction of more specific observables. 


5.2 The differential cross section e e + —> Z —> // 


The differential cross section for e _ e + —> /(si)/(s 2 ) with polarized final 
fermions can be written as a sum of spin independent, single-spin dependent 
and spin-spin correlation contributions as 


dcr 

dfy 


( s l> S 2) 


d<r° 

dff/ 


do - 5 

dsy (Sl ® S2) + 


dg c 

dfy 


(si <8> s 2 ) , 


(52) 


where Si (s 2 ) are the polarization vectors of the fermion (antifermion). We 
choose a reference frame in which the 2 -axis points in the direction of the 
outgoing fermion and the x—z plane contains both the incoming electron and 
the outgoing fermion. The polar angle Of is the angle between e~ and the out¬ 
going fermion. The azimuthal angle of the fermion can be integrated over, as 
the electron and positron are assumed unpolarized. The masses of the electron 
and positron will be neglected, as well as the quadratic terms in the anomalous 
couplings. The spin-correlation terms do not contribute when the polarization 
of only one of the two fermions is analyzed. 


In terms of the final fermion velocity /?, the dilatation factor 7 = Mz/Zmf = 
( 1 —/3 2 ) -1//2 in the overall c.m.s. and the polarization vectors 2 of fermion and 
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antifermion in their rest frames, the differential cross section e e + —> Z —> ff 
(exactly on the Z resonance) reads 


da° 

dcos Of 


a 2 /3N ( 


c 


x 


(v 2 e +ol)[vj(2-(3 2 shrO f ) 

+a 2 /3 2 ( 1 + cos 2 Of)} + 8 v e a e VfdfP cos Of 


(4 s w cw) 4 r| 

—8svt/CmRe(aJ) [(n 2 + a 2 )vf + 2 v e a e a,f/3 cos Of] 


(53) 


da 5 a 2 /3Nc 2s\yCw 

dcos Of (dswcw) 4 r| 

x [ ( s l + s 2)a;^+ + ( s l + ^ 2 ) 2 /^+ + ( S 1 + S^Z + 

+ (st - s;) x X_ + (si - s;) y F_ + (si - s * 2 ) z Z_ } , (54) 

with 


X + — l - [2v e a e v 2 + (v 2 + a 2 e )vfaf/3 cos Of 


{ 1S W C W 

+ 27 Re(aJ) [2n e a e n/(2 — /3 2 ) + {y 2 + a 2 )a//5cos6*/]| sin#/ , (55) 

Y + = — 2y/3lm(df) [2v e a e a,f + (v 2 + a 2 )vf(5 cos#/] sin Of , (56) 

Z + = —[(n 2 + al)vfdfP( 1 + cos 2 Of) + 2 v e a e (v 2 f + a,//3 2 ) cos#/] 

• s 'U' r U' 

+2Re(aJ) [4n e a e n/ cos #/ + (n 2 + a 2 )a//3(l + cos 2 0/)] , (57) 

X_ = 27 / 3 Im(dJ) [(n 2 + a 2 )vf cos Of + 2 v e a e a,ff3\ sin#/ , (58) 

Y_ = 2 'y/ 3 Re(ctf) [2 v e a e Vf + (n 2 + a 2 )a//5cos0/] sin#/ , (59) 

= — 2 f 3 Im(df) (n 2 + a 2 )vf sin 2 Of , (60) 

and 


d<r c 
dcos Of 


a 2 (3N } c 1 


SxxCxX d~ SyyCyy 4 “ S ZZ C Z 


(4 s w c w ) 4 r| 1 1 yy 

+( s xy + s yx)C xy + (-S^ + S 2a; + ( s yz + s zy)Cy Z 

+( S I x s V)xC yz + (si x s * 2 ) y C xz + (si x s * 2 ) z C xy } , 


with 


( 61 ) 


C xx = {yl + al)[vj (2 - ^ 2 ) - a^ 2 ] sin2 $/ 



— 8styCjyRe(a“) (u 2 + a 2 e )vf sin 2 Of , (62) 

Cyy = ~{v 2 e + a 2 e ){y) - a))f3 2 sin 2 9 f , (63) 

C zz = (v 2 e + a 2 )[u 2 (2 + /3 2 sin 2 9 f ) - a 2 f /3 2 ( 1 + cos 2 9 f )] 

+8 v e a e Vfdfl3 cos 9f 

— 2s\yCwR e ( a> f) [4(u 2 + a 2 )vf cos 2 9f + 8v e a e af/3cos9f} , (64) 

C+ y = 4s w c w /?Im(aJ) (u 2 + a 2 )a/ sin 2 9f , (65) 


= I ^[( v e + a l) v f cos0 f + 2u e a e u / a / /3] 

—4s^c^7Re(aJ) [{v 2 + a 2 )v/(2 — /3 2 ) cos#/ 


+2u e a e a//3] j sind/, (66) 

C+ = 4s w c W 7/?Im(aj) [(v 2 + a 2 )a/cos0/ + 2v e a e VfP] sin , (67) 

= 4s M /c M //5Re(ciJ) (u 2 + a 2 )v f sin 2 0/ , (68) 

C~ z = 4swCw'yPIm{df) [(v 2 + a 2 e )af/3cos9f + 2v e a e Vf]sin9f , (69) 

C~ z — — 4:SwCwl/3Re(ctf) [{y 2 e + a 2 )vf cos9f + 2v e a e af/3] sin#/ , (70) 


where s^- = Si ( s 2/ , i,j = x, y, z and AQC = 1 (3) for leptons (quarks). 

Notice that, from the SM at tree level, the final fermions are longitudinally 
polarized according to Eq. (57). There is also a SM contribution to the trans¬ 
verse polarization (within the collision plane) but it is helieity suppressed 
(see Eq. 55)). The transverse and the normal (to the collision plane) single¬ 
fermion polarizations are especially sensitive to the WDMs: the real part of 
the AWMDM (WEDM) contributes to the transverse (normal) single-fermion 
polarization (Eqs. (55) and (59), respectively), whereas the imaginary part of 
the AWMDM (WEDM) contributes to the normal (transverse) single-fermion 
polarization (Eqs. (56) and (58), respectively). The comparison of the polar¬ 
ization of both / and / isolates the CP-violating effects: the real part of the 
WEDM induces a difference in the / and / polarizations (Eq. (59)) orthogonal 
to the scattering plane whereas the imaginary part leads to a difference in the 
/ and / transverse and longitudinal polarizations (Eqs. (58) and (60), respec¬ 
tively). The spin-spin correlation terms also provide specific information on 
the dipoles according to their symmetry properties. 

Table 2 shows how the different terms of the effective vertex Zff (identified 
by their couplings) transform under C, P and T (y and a stand for whatever v t 
and cq, i — e, /).The properties of the transverse (sj x ), normal (s? ) and longi¬ 
tudinal (sj z ) polarizations are shown in Table 3. The cross section transforms 
as a real number under discrete transformations (it is C-, P- and T-even) 
and thus may only contain certain combinations of couplings and polariza¬ 
tions (Table 4). All these symmetries can be easily checked in the differential 
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Table 2 

Transformation properties of the different terms in the Zff effective vertex. 


Couplings 

P 

CP 

T 

V 

+ 

+ 

+ 

a 

- 

+ 

+ 

R e(aj) 

+ 

+ 

+ 

lm(aj) 

+ 

+ 

- 

R e(dj) 

- 

- 

- 

Im (dj) 

- 

- 

+ 


Table 3 

Transformation properties of single polarizations and spin-spin correlations. The 
axis are given by the direction of the fermion and the scattering plane as shown in 
the text. 


Correlations 

p 

CP 

T 

Sxxi Syy> Szzi 

[$xz “1“ $zx) 

+ 

+ 

+ 

($xy H - Syx)') 

( s yz + s zy) 

- 

+ 

- 

(si X S 2 )x,z 

- 

- 

- 

(si X S 2 )j, 

+ 

- 

+ 


Polarizations 

P 

CP 

T 

( S 1 + s 2)x,z 

- 

+ 

+ 

( S 1 + S2)y 

+ 

+ 

- 

(Si - S 2 )y 

+ 

- 

- 

(si — s 2 )a c,z 

- 

- 

+ 


cross section. 


5.3 Polarization analysis 


From the previous Section one concludes that the analysis of the spin-density 
matrix of the produced fermion pair (both the single-fermion polarization and 
the spin-spin correlation terms) is very important to disentangle the different 
anomalous contributions using appropriate observables. The polarization of a 
particle can be measured from the angular distribution of its parity-violating 
weak decay products. 

In the case of colored particles the QCD interactions play a crucial role. One 
needs that the quark decays before it can form hadronic bound states to 
analyze its polarization. This is not possible for light quarks and cannot be 
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Table 4 

Allowed combinations of couplings and polarizations in the cross section. 



SM 

Re(aJ) 

Im(ay) 

R e(dj) 

Im (d,J) 

no spins 

v 4 , v 2 a 2 , a 4 

S 2 

v , va 

— 

— 

— 

x + , z + 

v 3 a, vo? 

v 2 a, a 3 

— 

— 

— 

T+ 

— 

— 

^ 9 

r, va 

— 

— 

X-, z_ 

— 

— 

— 

— 

S 2 

v , va 

y_ 

— 

— 

— 

v 2 a, a 3 

— 

C'xx ? C 'yyi C'zz ? ^xz 

v 4 , v 2 a 2 , a 4 

^ 9 

v°, va 

— 

— 

— 

w xy » ^yz 

— 

— 

v 2 a, a 3 

— 

— 

n- n- 

KJ xy'< yz 

— 

— 

— 

Q 9 

v' , va 

— 

c- 

— 

— 

— 

— 

v 2 a , a 3 


done with reliable precision for the h quark [31]. The angular distribution of 
the jets in the 2-jet exclusive decays of the Z does not carry any information 
on chirality-conserving effective couplings [48] and cannot probe CP, T or 
CPT invariance even in the case when the flavor of the jets can be tagged 
[49]. CP-violating effects can be studied in 2-jet inclusive and 3-jet exclusive 
decays of the Z [50-52], The t quark is the only one heavy enough to decay 
before hadronization but it deserves a different treatment as it cannot be 
pair-produced at the Z peak. 

The r lepton case will be studied below as a prototype for polarization anal¬ 
yses. In the case of r pair production the expressions are simpler [27,29], as 

X 1 

v = v e = v T = -- + 2, a = a e = a T = -- . (71) 

For the single polarization analysis below the direction of flight of the r must 
be reconstructed. It has been shown [53] that this is possible if both r 1 * 1 decay 
semileptonically (at least one hadron in each decay) and the energies and 
tracks of both hadrons are determined (using microvertex detectors). Consider 
the r-pair production and their ulterior decay, in the case where at least one 
of them decays semileptonically (r“ —» hfv y or r + —> fbfvr)- The direction of 
emission of the final charged hadrons in the semilcptonic decays acts as a r ± 
spin analyzer. Using the narrow width approximation [54], one can write the 
cross sections for the whole process, in terms of the polar and azimuthal angles 
of the charged hadrons in the t± rest frames in a very easy way: substitute 
the polarization vectors s± of the t 31 by =F«/i± • 91h±, respectively, where is 
the unit vector of the corresponding charged hadron in the r ± rest frame. The 
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parameter ah is a measure of the spin-analyzing power of the corresponding 
decay. It is maximal, a n — 1, for the decay r —> ttu t and it is given by 
a p = ( m 2 — 2m 2 )/(m 2 + 2m 2 p ) ~ 0.46 for r —> pu r . 

Integrating over the angular distribution of the decay products of the r whose 
polarization is not analyzed, the spin-spin correlations between the two decay 
branches vanish. In this case, one obtains: 


dcr 


dn T -dn h T 47 t 


= —BR(r T -> h T ->■ hfX) 


x 


do -0 


± 


a 2 /3 2 s\yCw 


dO T - (4 swcw) 4 r| ^ 

+ (y+ ± y_) sin 6 h T sin (j) hT + (Z + ± ZJ) cos 9 h ^ 


(X + ± XJ) sin 9 h ^ cos (f)^ 

(72) 


This expression can be further simplihed integrating over the polar angle of 
the analyzing hadron, still keeping valuable information on all the anomalous 
couplings [27], 


da 


d cos 0 T -d(f>hi 


= BR(r F -> i' r ) )BR(T ± -> h^X) 


do - 0 7 Ta 2 f3 2 s w c w 

dfl T - (4 s w c w ) 4 T| ^ 

+ (y+ ± Y_) sin 6 h ^ sin (j) hT 


Oik =F ■ { {x + ± XJ) COS (j) h T 

(73) 


In the single-polarization distributions above, the CP-even and CP-odd con¬ 
tributions cannot be disentangled, as the CP-conserving (CP-violating) ef¬ 
fects represented by the real part of the AWMDM (WEDM) may be faked 
by absorptive effects (unitary corrections due to the exchange of virtual par¬ 
ticles). For instance, as mentioned above, the transverse polarization reads 
both Re(aj) and Irn(riy) and the normal polarization reads both Re(dJ) and 
Im(ay). In principle, one needs to compare the polarizations of both pair- 
produced fermions to eliminate the component of unwanted CP parity. Nev¬ 
ertheless, exploiting the (accidental) large suppression of the vector coupling 
of charged leptons one can still define a set of observables based on single 
polarizations which are sensitive to the dispersive and absorptive parts of the 
t weak dipole moments [27,29,30]. A genuine CP-sensitive observable must 
involve the two pair-produced fermions. This can be done comparing the two 
decays (at the price of half the statistics) or using spin-spin correlations, which 
manifest in angular correlations among their decay products. 
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5.3.1 Observables from single r polarization 

We list below the observables sensitive to the WDMs proposed by Bernabeu 
et al. [27,29,30] based on single r polarization.^] They consist of azimuthal 
asymmetries on the analyzing hadron direction and demand the reconstruc¬ 
tion of the r frame. They exploit the helicity flipping character of the dipole 
moments and the symmetry properties of the anomalous terms. The present 
limits on the r AWMDM, that will be quoted later, are based on a couple of 
these asymmetries. 


Real part of the AWMDM Using the distributions (73) one select events 
sensitive to the transverse polarization through 


cr ^(+) 



3tt/2 


J d cos 9 r - J d(j) h =F + j d cos 0 T - J d 

0 —7t/2 — 1 7t/2 

1 37t/ 2 o 7r/2 

J dcos#,-- J dcfhT + J d cos 9 r - J dcf)^ 


-71-/2 


dcr 

dcos 9 T -d(j)f l T 
da 

dcos 9 r -dfhT ’ 

(74) 


from which one can construct asymmetries especially sensitive to the dispersive 
part of the AWMDM 


_ ^c(+) ~ 
““*£(+) + *£ (-) 


a h T 


a 


-s w cw 


v 2 + a- 


-+ 2 7 Re«) ± 2 7 -Im« 

7 SVV C W a 


(75) 


The tree level contribution is helicity suppressed. A possible absorptive con¬ 
tribution from the WEDM is suppressed by v/a and can be ignored. Notice 
that it can otherwise be eliminated using a genuine CP-even observable that 
compares both t ± decaying into the same kind of hadrons (a = a h + = «/,-): 
7f C c = \{A~ C - A+) oc Re«). 


Imaginary part of the AWMDM The following asymmetries related to 
the normal polarizations are sensitive to the absorptive part of the AWMDM: 


13 Here we present exact expressions in agreement with [27,29,30] only for v <C a. 
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A possible contribution of the WEDM is suppressed by v/a. (Again a true 
CP-even observable like A s = |(A“ — A+) is completely free from the WEDM 
dependence.) 


Real part of the WEDM From events sensitive to the normal polarization 
of the r’s 


<(+) = 



r i 

J d cos 9 t - 

- o 
r 1 

J dcos 0 T - 

o 


0 

2 IT 


u 

d0/j T + J dcos 0 7 

-l 

o 

d 4>hT + J dcos 8 r 



da 

dcos 0 T -d(f) h T 
da 

dcos 6h-d0/jT 


(77) 

(78) 


the following asymmetry isolates the real part of the WEDM with a v/a 
suppressed influence from the absorptive part of the AWMDM: 


= “ ah ^ SwCw lF^ 


Re(d/ 


=F -Im(a* 
a 


• (79) 


The true CP-odd observable Acp = |(A SC + A+) is completely free of the 
contribution from AWMDM, as expected. 


Imaginary part of the WEDM A genuine CPT-odd and CP-odd observ¬ 
able based on single r polarizations could be constructed from (75) although 
it is numerically very suppressed: 


A 


c p _ 

CC 



^ h 

— 2 " 7 s wCw 


v 2 + a- 


-Im 


(80) 
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5.3.2 Observables from spins of both fermions 


In this Section, the polarizations of the two fermions are analyzed simulta¬ 
neously. Usually this kind of observables involve spin-spin correlations (i.e. 
angular correlations of the final particles in the decays) although in some case 
only single fermion polarizations might contribute. We restrict ourselves to 
observables sensitive to CP violation. As already mentioned, only when both 
fermions (particle and antiparticle) are considered one can get genuine indi¬ 
cations of CP violation. The methods described below are also suitable (and 
extensively applied) to search for CP violation in e + e _ —> // away from the 
Z resonance, where other CP-violating form factors, besides d z (s), can be 
relevant. Our starting point is an initial CP-even eigenstate, in the c.rn.s., as 
e + and e~ beams are assumed unpolarized. Otherwise a non-zero expectation 
value of a CP-odd observable would not imply, in principle, CP violation. This 
prescription can be relaxed in some situations: 

• For the case where the e + and e~ are transversely polarized, opposite in 
direction but equal in magnitude, the initial state is not a CP eigenstate 
but certain CP-odd observables still provide a signal of CP violation [55]. 

• Recently, also the possibility of using longitudinally polarized beams has 
been proposed. Again in this case the initial state is not a CP eigenstate and 
therefore a CP-odd and CPT-even observable (sensitive to Re(cTJ’ z )) can 
be contaminated by CP-conserving, absorptive (T-odd) effects; similarly a 
CP-odd, CPT-odd observable (sensitive to lm(dy z )) can be contaminated 
by CP-conserving, dispersive (T-even) effects. Nevertheless it has been ar¬ 
gued that these contaminations are at the per mill levelpH [56,58] and hence 
far from the sensitivity of a high luminosity e + e _ collider: true CP-violating 
effects, if any, would not be faked. 

Consider e + (p + ) + e“(p_) —> r + (k + ,s + ) + r _ (k_,s_) with t~ —> a(q_) + 
neutrals and r + —» 6(q + ) + neutrals. Then the momenta and polarization 
vectors in the overall c.m.s. transform as follows: □ 


14 In the SM, at leading order in the electroweak couplings, the fermion pair produc¬ 
tion proceeds through an intermediate CP-even eigenstate (7 or Z) and hence, to 
one-loop, only electroweak box contributions would be relevant for the absorptive 
effects, but they are estimated to be negligible [56]. For the dispersive CP-conserving 
contributions, they can only arise in the SM from bremsstrahlung off the initial e~ 
and e + , but this effect is also small [57,56]. 

15 These vectors refer to a fixed frame. The transformation properties of s± are not 
in contradiction with Table 3 as s]j_ refer to the frame given by the scattering plane, 
with axes: z = k^, y = k_ x p_/|k_ x p_| and x = y x z. The symmetries of the 
different components of s]j_ are thus connected with the ones of k_, p_ and s±. 
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( 81 ) 


CP: p± 

“Pt 

= p± 

CPT : p± 

k± 

H- 

\ 

T 

= k± 

k± 

s± 



s± 

q± 

-»■ -q=F 


q± 


Pt 

k T 

— S =F 

q F 


-p± 
—k± 


From the unit momentum of one of the r in the c.m.s. (e.g. k + ) and both r 
polarizations (s-t) a basis of linearly independent spin CP-odd tensor observ¬ 
ables can be constructed [59]. They are classified according to their “0 =CPT 
parity” given by (—l) ra *4. e = p^A with n being the rank of the observable. 
We show here, for later reference, the scalar observables: 


A {1) = k + • (s+ - s_) [CPT-odd] (82) 

A^ = k + ■ (s + x s_) [CPT-even] (83) 

and two vector observables: 

^ (3) = s + - s_ [CPT-odd] (84) 

A^ = s + x s_ [CPT-even] (85) 

Other spin observables are based on them. For instance [60] 

(p+ x k + ) • (s+ — s_) [CPT-even] ( 86 ) 

is just a projection of A <:i> sensitive to the CP-odd combination of the trans¬ 
verse polarizations, (s_ — s+) y . Similarly, A (l) is nothing but (s_ — s + ) z (see 
Table 3). 


The spin observables are related to more realistic (directly measurable) mo¬ 
mentum observables based on the momenta of the r —> 1 charged prong decays 
[61]. The dimensionless ones are easier to measure, for instance the scalar ob¬ 
servables^] [62,63]: 


ii = p+ • ,^ + X [CPT-even] (87) 

|q+xq_| 

A 2 = p+ ■ (q+ + q_) [CPT-odd] ( 88 ) 

or the CP-odd traceless tensors [59,62,63]: 

16 Notice that A\ is related to W 2 ) and A^\ The observable A 2 is related to *4®. 
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Tij - (q+ q-)i 1 \ J + (*<">■ j) 

[CPT-even] 

(89) 

q+ x q-l 



Qij = (q+ + q-)»(q+ - q-)j + (*<-»• j) 

[CPT-odd] 

(90) 

The usual decay channels are r —> £ugu T (£ = 

e, p), t —> 7iu T and r - 

-+ pv T 


ttti 0 v T (they amount to about 70% of the r decays). Notice that the recon¬ 
struction of the r frame is not necessary for the momentum observables and 
hence the restriction to semileptonic r decays is no longer mandatory. As a 
final comment, the observables and Q tJ do not involve angular correlations 
as they could be measured considering separate samples of events in the re¬ 
actions e + e _ —> aX and e + e~ —> aX. Nevertheless it is convenient to treat 
them in an event-by-event basis [62], 

One may obtain additional CP-odd observables from combinations of the 
ones given above by multiplying them with CP-even scalar weight functions 
to maximize the sensitivity to CP-violating effects (optimal observables). Ne¬ 
glecting quartic terms in the CP-violating couplings, the differential cross 
section can be written in terms of the the real and imaginary parts of the 
EDM and WEDM (A*, i — 1, ...,4), da = dao + A,da}. It has been shown 
[64,65] that the observables given by O, = da^/dao have maximal sensitivity 
to the CP-violating parameters A*. 


5-4 Experimental limits from these observables and comparison with theoret¬ 
ical predictions 


One can give a rough estimate of the sensitivity of a certain experiment to 
put limits to the form factors on a simple statistical basis. This simple ap¬ 
proach provides always optimistic results. More refined analyses require the 
implementation of the detector performance. 

Let A be an asymmetry proportional to an anomalous form factor F, A — c-F. 
The accuracy to determine this asymmetry from a sample of N events in an 
experiment is given by A(M) = n\J (1 — A 2 )/N ~ n/yfN, as the asymmetry 
is expected to be very small, n is the number of standard deviations (s.d.) 
related to the confidence level of the measurement. If no significant effect is 
observed in the experiment the accuracy translates into an upper bound given 
by |F| < n/(\c\VN). In our case, N = (# of Z)-BR(Z -> t~t + )-BR(t~ -> 
A)-BR(r + —> B) if the decay channels A and B are considered as spin ana¬ 
lyzers. 

The expectation value of a CP-odd observable O is given byf^] (0) o6 = 
\((0)ab + (0)ba) = C ab ■ with (°)ab = $ ai®/1 ^ab where the T decay 

17 This relation holds only at the Z peak as the on-shell amplitudes contain only 
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channels a and b are selected. Of course any phase space cut must be CP- 
symmetric. The accuracy of a measurement based on a sample of N events is 
A (O)ab = nsj (( 0 2 ) ab - (0) 2 ab )/N ~ nsj ( 0 2 ) ab )/N , as the CP-odd effects are 
expected to be very small. The differential cross section is d a a i ~ da®^+dcdy p , 
where the second term is CP-odd (linear in d™). In practice, the average ( 0 2 ) a b 
is given by the SM. If the measured value is (compatible with) zero, an upper 
limit given by d™ < nsj(0' 2 ) a b) /(| c 0 j } \ \J~N) with the confidence level corre¬ 
sponding to n s.d. . 

It is customary to dehne the ratio r = (O)/\J(O 2 ), which is a measure of 
the sensitivity of the CP-odd observable [56,58,60]. The optimal observables 
are intended to maximize |r|. Given a sample of IV = 2N + = 2AL events, for 
the process and its CP image, the corresponding signal-to-noise ratio is then 
S = y/N\r |. This is nothing but the statistical significance of the signal: S — n 
corresponds to an n s.d. effect. 

The present experimental limits are based on measurements compatible with 
zero for the AWMDM and WEDM: 

• Recently a measurement by the L3 detector at LEP of the weak magnetic 
dipole moment of the r lepton has been performed, for the first time, using 
the azimuthal asymmetries for the channels r —>• nu T and r —> pu T [66]. Also 
a limit to the weak-electric dipole moment was reported. The results are 

|Re«)| <6 x 1(T 3 , (91) 

|Im«)| <2 x 1(T 2 , (92) 

|Re«)| <4 x 10~ 17 ecm = 7.2 x 10“ 3 p T (93) 

(95 % C.L.), an order of magnitude less stringent than the expectations in 
[30], where all the semilcptonic decay channels (52% of the total decay rate) 
are considered. The d™ result is not competitive with that obtained by the 
CP-odd tensor methods, given below. 

• For the WEDM the present limits were obtained by the OPAL collaboration 
using optimal observables [67] with 95% C.L.: 

|Re«)| < 5.6 x 10' 18 ecm = 10~ 3 /j r , (94) 

|Im«)| < 1.5 x 10- 17 ecm = 2.7 x 10“ 3 p T . (95) 

These experimental limits are not in conflict with the theoretical expectations 
of the SM and the MSSM. They both predict a real (imaginary) part for the r 
AWMDM at least three (four) orders of magnitude smaller (19, 32) than the 
present experimental sensitivity (91, 92). Closer are the r WEDM expectations 
in the MSSM: the real part can reach values of the order of 3 x 10” 20 ecm 

one CP-odd form factor, the WEDM. 
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in the large tan/3 scenario (39) roughly two orders of magnitude below the 
experimental reach (94,95). 


6 Observables off the Z resonance: the t quark case 


One needs to go beyond the Z peak to produce t quark pairs and this implies 
to account for several new features: 

(1) The photon-exchange diagram is no longer suppressed in this regime and 
therefore one needs to separate the contributions of electromagnetic and 
weak form factors. It is then enough to employ as many independent 
observables as necessary. For instance, looking for CP-violating effects, 
one may use T tJ and A\ to isolate the real parts and Q l3 and A 2 to get the 
imaginary parts of dj(s) and df(s) (cf. Ref. [63] for the r and Ref. [62] 
for the t quark). 

(2) Not only vertex- but also box-diagrams correct the tree level process 
to one loop. The former corrections are parameterized by the electro¬ 
magnetic and the weak vertex form factors but the latter demand the 
introduction of new form factors according to the more general topology 
of the process. 

(3) In general, the dipole form factors a j (s) and dj (s) are not gauge invariant 
to one loop for an off-shell V. Exceptions are the CP-violating ones 
because they do not receive contributions from vertex corrections in which 
gauge or Goldstone bosons are involved. This at least guarantees the 
gauge-parameter independence for the class of R% gauges. 

In any case, for a realistic theory, one expects that any CP-odd observable 
will depend not only on the CP-violating effects due to vertex corrections 
parameterized by the dipole form factors but also on possible CP-violating 
box contributions. We concentrate on supersymmetric CP violating effects in 
f-pair production at e + e _ colliders. CP conserving MSSM one-loop contribu¬ 
tions to e + e _ —> ff are discussed in [ 68 ]. To our knowledge, only the electric 
and the weak-electric form factors have been considered so far to parameterize 
these effects. □ Our purpose is: 

• To briefly review a specific set of CP-odd observables in f-pair production. 

• To evaluate the expectation value of several observables in the context of 
the MSSM with complex parameters to one loop. 

• To compare the contribution from the electric and weak-electric form fac¬ 
tors to these observables with the CP-violation effects coming from the 
box contributions. In this way we illustrate, with the MSSM as a concrete 

18 A similar analysis for hadron colliders has been recently presented in Ref. [69]. 
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example, the importance of the box effects. 


6.1 CP-odd, observables 


The spin effects can be analyzed through the angular correlation of the weak 
decay products, both in the nonleptonic and in the semileptonic channels 


t{ k+) -> &(qb)^h ad (qx) , (96) 

*(k+) -> b(q h )£ + (q + )u e (£ = e,n,T) (97) 


and the charged conjugated ones. □ As m t > M w + nib, the 1 quark decays 
proceed predominantly through Wb. Within the SM the angular distribution 
of the charged lepton is a much better spin analyzer of the t quark than that of 
the b quark or the W boson arising from semileptonic or nonleptonic t decays 
[71]. Usual CP-odd observables are: 


(1) The analogous to the ones for the t + t~ case [62] (Eqs. (82-90)), including 
the optimal observables [64,72], obtained from their multiplication with 
appropriate CP-even scalar weight functions depending on the energies 
or momenta of the final state particles that maximize the sensitivity to 
CP-violating effects. 

( 2 ) As the t quark is a heavy fermion, the CP conjugate modes tidL and t R t R 
are produced with a sizeable rate. This allows to construct the following 
CP-odd asymmetry [73,30] 


#QUU.) ~ #(t R t R ) 
#(t L i L ) + #(t R t R ) 


[CPT-odd] 


(98) 


This asymmetry is related to the one that can be measured through the 
energy spectra of prompt leptons in the decays t —>■ W + b —> l + V(b and 
its conjugate. The W + is predominantly longitudinally polarized and, 
assuming the standard V—A interaction, the b quark is preferentially left- 
handed. The W + is mostly collinear with the t polarization and so is the 
£ + anti-lepton. Above the tt threshold a £ + coming from a t R has more 
energy, due to the Lorentz boost, than one produced in a U decay. The 
same happens for the conjugate channel, the l~ from a t l is in average 
more energetic than the one from the t R . Therefore, in the decay of the 
t R t R the anti-lepton from t R has a higher energy E + , while in the decay 
of the pair tide the lepton from tj, has a higher energy E_. Thus the 
asymmetry A is sensitive to the energy asymmetry of the leptons [49,74], 


19 The leading QCD corrections to e + e —* tt, that include a gluon emission, have 
a very small effect on the t spin orientation [70]. 
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Fig. 13. Generic box graphs for the MSSM at one loop. 


6(E + -EJ)~ 6(E _ - E+) [CPT-odd] (99) 

where 9(x) is the step function. Using the b quark as spin analyzer, a 
similar asymmetry based on b and b energies has also been proposed [75]: 


A£ p = Af - A? [CPT-odd] 

,E — W(6) > 3>) ~ < -E’o) 

W11 m = #7iw7 > -Bo) + si < -Bo) 


( 100 ) 


where E 0 = s/sfrn^ — M^) / Arri^ is the average energy of the b or b quarks. 
(3) The T-odd triple product correlations [62,60,56,58,36]: 


A^ = (T + )-(T_) [CPT-even] 

with T + = (qiq 2 q 3 ) = (qi X q 2 • q 3 ) and 71 = CP(T + ) (101) 


and T-even correlations [60,75]: 

= (V+) - (X>_) [CPT-odd] 

with V + = qi • q 2 and = CP(77 + ) (102) 


and their dimensionless versions, where qi 23 can be any of the 3-momenta 
in e + e _ —> tt or the t (t) decay products and A 2 are just two of them). 
These asymmetries have been proven suitable to investigate whether the 
source of CP violation resides in the tt production or in their decays. In 
particular, the ones based on (q+qxP+) and qx • q+ are sensitive to CP 
violation in tt production while the asymmetry from (q+qt,P+) (close to 
threshold) traces CP violation in semilcptonic t and t decays [60]. 

We will avoid the possible CP violation in the t or t decays by considering 
(ideal) spin observables. In this way we face the evaluation of the influence on 
the CP-odd observables of the vertex and box diagrams without any interfer¬ 
ence of other CP-violating effects in the t decays. The expectation value of 
some of the realistic momentum observables given above will also be presented 
(assuming SM top decays) for comparison with experimental capabilities. 
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Table 5 

CP-odd spin observables and the coefficients for the expectation value of their inte¬ 
grated version at ffis = 500 Ge V, where only the CP-violating dipole form factors 
are taken into account. 


i 

CPT 

Oi 

a 

b 

Cl 

C2 

C3 

c 4 

1 

even 

(si - s^ 

N| 

N| 

0.526 

0 

1.517 

0 

2 

even 

(si X sD* 

NT 

LT 

-0.465 

0 

-0.061 

0 

3 

even 

(si X s* 2 ) z 

NT 

TT 

0.708 

0 

0.144 

0 

4 

odd 

(si - S 2 ):r 

TT 

Ti 

0 

0.930 

0 

0.123 

5 

odd 

(si - S* 2 ) z 

LT 

Li 

0 

-1.417 

0 

-0.287 

6 

odd 

(S*l X S 2 ) y 

LT 

Ti 

0 

0.263 

0 

0.758 


6.2 MSSM full contribution to CP-odd observables 


The one-loop differential cross section for polarized i-pair production in the 
MSSM [68] involves the box diagrams indicated in Fig. 13 besides the vertex 
graphs of Fig. 1. The class with vector boson exchange contains only SM con¬ 
tributions ([eZZt], [u e WWb]) whereas the other one is purely supersymmetric 
([ex°x°t], [H'eX ± X ± b]) and provides CP violating effects. The box diagrams 
with Higgs boson exchange are proportional to the electron mass and can be 
neglected in the whole calculation (they are CP even in any case). When we 
refer to one-loop calculation in the following, the QED as well as the standard 
QCD corrections to the tree level process are excluded: they need real photonic 
and gluonic corrections to render an infrared finite result and constitute an 
unnecessary complication as they are CP-even and do not affect qualitatively 
our conclusions. 


6.2.1 Spin Observables 

Consider the process e + (p)+e _ (—p) —> t( k, Si)+t( —k, s 2 ) (pair-production of 
polarized t quarks). A list of CP-odd spin observable^ 7 ] classified according 
to their CPT properties is given in Table 5. Their expectation values as a 
function of s and the scattering angle of the t quark in the overall c.m. frame 
are given by 


20 Notice that Ci, 3,5 are equivalent to (86,82, 83) respectively and the rest are some 
other projections of (84) and (85). s| 2 are the t and t polarizations defined in their 
rest frames. 
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(O) ab = 


2d a 


E + E 

|_s£,S2=±a,±b s i> s 2 == tb,±aJ 


dcr(s*,s;) O , 


da= ]T d(J ( s l,4)- 

d=Si ,d=S2 


(103) 

(104) 


The directions of polarization of t and i (a, b) are taken normal to the scatter¬ 
ing plane (N), transversal (T) or longitudinal (L). They can be either parallel 
(I) or antiparallel (j) to the axes defined by z = k, y = k x p/|k x p| and 
x = y x z. Notice that 


(O 2 ) ab = ^- E dcr(si,s;) (9 2 (105) 

a(J sj,sj=±a,±b 

and the same for the average of any CP-even quantity. If the information of 
the t scattering angle is not available one can consider integrated observables 


< 0 >ab 


1 

2a 


E 

sps*=± a ,±b 


E 


spS2==tb,±a 


a s 


1) a 2/ 


O . 


(106) 


The contributions to the CP-odd observables are linear in the t EDFF and 
WEDFF and in the CP-violating parts of the one-loop box graphs. The shape 
of the different dipole contributions to these observables is depicted in Fig. 14. 
The coefficients of the dipoles in the linear expansion of the integrated spin 
observables, 


(O)ab = (ci Re[d?(s)] + c 2 Im[d7(s)] + c 3 Re[df («)] + c 4 Im[df (s)]) 

(107) 

are shown in Table 5 for their integrated version at y/s = 500 GeV. They are 
model independent. 

In Fig. 15 we compare the contributions of dipoles and boxes to the spin 
observables, for two reference sets of parameters: the first was given in (48) 
and is the one that has been shown to enhance the dipole effects; the second 
is given by 

Reference Set #2 : tan/3 = 1.6 

M 2 = H = rriq = \m t LR \ = \m h LR \ = 200 GeV 
W = <Pi = l Pb = E 2 1 (!08) 
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( 1 , 0 , 0 , 0 ) 


( 0 , 1 , 0 , 0 ) 



cos 0 


cos 0 


Fig. 14. Dipole contributions to the expectation values of the spin observables at 
= 500 GeV for different values of (Re[df], Im[dJ], Re[df], Im[dfj) in /at units. 


where, compared to (48), only the CP-violating phases have been changed. 
The shape of the solid and dashed curves is the same in all cases, as expected. 
Their different size is due to the contributions to the normalization factors 
coming from self-energies, A(W)MFFs and other CP-even corrections. The 
plots show that for the Set ff 1 the MSSM box graphs contribute in general 
to CP violation in the process e + e~ —> tt by roughly the same amount and 
with a different profile than the EDFFs and WEDFF of the MSSM. For nearly 
all cases this eventually results in a large reduction of the value of any CP- 
odd observable with respect to the expectations from the dipoles alone, as it 
will be shown below. Other sets of SUSY parameters, that do not enhance 
the dipole contributions, can provide instead smaller dipole form factors but 
larger observable effects. This is the case of the Set #2, that yields much 
smaller values for the real parts of the dipole form factors at a/s = 500 GeV 
than we had before in Eq. (50), 


df = (0.547 - 1.889 i) x 10~ 3 /x u 


(109) 
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Reference Set #1 



-1 -0.5 0 0.5 

cos 8 


-1 -0.5 0 0.5 

cos 6 


-1 -0.5 0 0.5 1 

cos 8 


Reference Set #2 



cos 8 


cos 8 


cos 8 


Fig. 15. Expectation value of the spin observables [in 10 -3 units] for the two ref¬ 
erence sets of SUSY parameters, assuming for the cross section: (i) tree level plus 
contribution from (W)EDFFs only (solid line); (ii) one loop including all the vertex 
corrections and the self-energies (dashed line); (in) complete one loop (dot-dashed 
line). 
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Table 6 

Ratio r = (0)/ \J ( 0 2 ) [in 10~ 3 units] of the integrated spin observables at -^/s = 500 
GeV for the two reference sets of SUSY parameters. The left column excludes the 
box corrections and the right one comes from the complete one-loop cross section 
for e + e~ —> tt. 


i 

CPT 

Oi 

a 

b 

Set ffl 

Set #2 

1 

even 

(si - S * 2 )y 

Nj 

N| 

1.216 

-0.231 

0.207 

-1.394 

2 

even 

(si X 

N| 

LT 

-0.755 

-0.489 

-0.053 

0.318 

3 

even 

(si X S * 2 )z 

N| 

TT 

1.184 

0.625 

0.090 

-0.598 

4 

odd 

(si - s* 2 ) x 

TT 

T| 

-1.230 

-1.421 

-1.888 

-2.217 

5 

odd 

(si - s^) 2 

LT 

Li 

2.550 

2.739 

3.823 

4.216 

6 

odd 

(Si X S l)y 

LT 

T| 

-1.683 

-1.751 

-2.050 

-2.216 


df = (0.362 — 2.319 i) x 10 -3 ^t- (110) 

As a consequence, the CPT-even observables receive their contributions mainly 
from the box diagrams (Fig. 15). 

The previous arguments are reflected in Table 6. There we show the ratio r = 
( O) /( O 2 ) for the integrated spin observables. We compare the result when 
only the self energies and vertex corrections are included (left column) with the 
complete one-loop calculation (right column). The shape of the observables 
as a function of the t polar angle (Fig. 15) is such that the signal for CP- 
violation r is still not very large for a couple of CPT-even observables with 
the Set f [2 but it is much better for all the others. This illustrates that the 
dipole form factors of the t quark are not sufficient to parameterize observable 
CP-violating effects and the predictions can be wrong by far. 

As final comments, notice that O 5 is nothing but the asymmetry of Eq. (98), 


iQ \ _ #(U^l) ~ #(t R t R ) 
#(Wl) + ) 


( 111 ) 


Attending to Table 5 it happens to be the most sensitive observable to the 
imaginary part of the EDFF and also the best one to test CP violation for our 
choice of SUSY parameters (Table 6). The observables (s£ — sJj) can still give 
information on CP violation when the polarization of one of the t quarks is not 
analyzed. □ The sensitivity of the single-spin polarization to CP-violation is 

21 Notice that anyway a comparison between two samples, one with polarized t and 
the other with polarized t, is necessary to build the genuine CP-odd observable of 
Eq. (104). 
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Table 7 

Ratio r [in 10 -3 units] of the momentum observables (87-90) and their ratio r at 
yfs = 500 GeV for three different channels: t + t—* a(q+) + c(q_) + X, given the 
Set #2 of SUSY parameters (108). The left column includes only the t (W)EDFF 
corrections and the right one comes from the complete one-loop cross section for 
e + e~ —> tt. 


CPT 

O 

b- 

b 

l- 

-b 

e- 

-t 

even 

M 

-0.036 

0.242 

0.030 

-0.202 

0.068 

-0.467 

odd 

A-2 

0.270 

0.304 

-0.180 

-0.204 

-0.501 

-0.812 

even 

T33 

-0.006 

0.042 

0.021 

-0.140 

-0.037 

0.248 

odd 

Q33 

0.486 

0.542 

-0.335 

-0.374 

-1.274 

-1.420 


indeed worse: ((si - s* 2 ) xy J) ah 


2((si) w )a for a(b)= NT(j),TT(j),LT(j). 


6.2.2 Momentum Observables 

Consider now the decay channels labeled by a and c acting as spin analyzers 
in t +1 —> a(q+) + c(q_) + X. The expectation value of a CP-odd observable 
is given by the average over the phase space of the final state particles, 

{0)ac — ~ [(C^)ac+ {0)cs\ = -- f [d(T a c + d <7 C a\ O , (112) 

^ ac ^ 

where both the process (ac) and its CP conjugate ( ca ) are included and 


Oac J d(J ac J d(J ca 5 (113) 

in full analogy with Eqs. (103,104). The differential cross section for f-pair 
production and decay is evaluated for every channel using the narrow width 
approximation. 

In Table 7 the ratio r is shown for three different decay channels at \fs = 500 
GeV and some realistic CP-odd observables involving the momenta of the 
decay products analyzing t and t polarizations in the e + e _ c.m.s. (laboratory 
frame). As expected the leptonic decay channels are the best t spin analyzers 
but the number of leptonic events is also smaller. The Reference Set ff 2 has 
been chosen. As expected, the dipole contributions (left columns) to the CPT 
even observables are very small for this choice of SUSY parameters but the 
full expectation values (right columns) are larger. 

Concerning the experimental reach of this analysis, the statistical significance 
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of the signal of CP violation is given by S = \r\\/N with N = eCa t iBK(t —> 
a)BR(£ —> c) where e is the detection efficiency and C the integrated luminosity 
of the collider. The branching ratios of the t decays are BR ~ 1 for the b 
channel and BR ~ 0.22 for the leptonic channels {£ = e, /i). At y/s = 500 
GeV the total cross section for t-pair production is ~ 0.5 pb and the NLC 
integrated luminosity C ~ 50 fb” 1 [21], Assuming a perfect detection efficiency, 
one gets y/N ~ 160, 75, 35 for the channels b — b, £ — b, £ — £, respectively. 
With these statistics, values of |r| ~ 10 -2 would be necessary to achieve a 1 
s.d. effect, which does not seem to be at hand in the context of the MSSM as 
Table 7 shows. Nevertheless, the ratio r should improve considerably (~30%) 
using optimal observables (suitable to enhance the dipole effects) and also 
with (longitudinally) polarized beams [56]. 


7 Summary and conclusions 


The one-loop expressions of the dipole form factors of fermions in terms of ar¬ 
bitrary complex couplings in a general renormalizable theory for the ’t Hooft- 
Feynman gauge have been given. The CP-violating (-conserving) dipole form 
factors depend explicitly on the imaginary (real) part of combinations of the 
couplings. 

As an application, the weak-magnetic and weak-electric dipole moments of the 
r lepton and the b quark (defined at s = M§) and the electric and the weak- 
electric dipole form factors of the t quark (for s > 4m 2 ) have been evaluated for 
the MSSM with preserved R-parity and non-universal soft-breaking terms. A 
version with complex couplings of the MSSM has been implemented in order to 
get one-loop CP-violating effects. The quark dipoles depend on three physical 
CP phases and the lepton ones only on two since there is no mixing in the 
scalar neutrino sector. The supersymmetric parameter space has been scanned 
in search for the largest contributions to the dipoles. Both the r and b dipole 
moments are enhanced for large tan f3. There the AWMDMs can be one order 
of magnitude larger than in the electroweak SM (in particular the real parts, 
(P(1CT 5 )) but still a factor five below the standard QCD contribution for the b 
case. The WEDMs can be instead twelve orders of magnitude larger than the 
tiny three-loop prediction by the SM and reach 0( 10” 20 ) ecm for the r and 
(P(1CT 19 ) ecm for the b. Conversely, the t dipole form factors are larger in the 
low tan/3 scenario. The values depend strongly on the interplay between the 
energy at which they are evaluated and the set of supersymmetric parameters 
used as inputs. The real and imaginary parts can reach values of similar size: 
at a/s = 500 GeV, the t EDFF and WEDFF are of (P(10~ 19 ) ecm. 

At the Z resonance the electromagnetic dipole form factors are irrelevant 
and the weak dipole form factors (weak dipole moments of the pair-produced 
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fermions) are gauge invariant objects that can also be related to physical 
observables. Such observables are based on the polarization analysis of the 
fermions, equivalently on the angular and energy distributions of their decay 
products. The cross section for polarized fermion production in e + e _ collisions 
at the Z peak has been given in terms of the AWMDM and WEDM, as well 
as a review on observables able to discriminate the dipole moments for the r 
case (the polarization analysis does not appear feasible for the b case). The 
present experimental limits based on them have been compared to the SM 
and the MSSM expectations. The experimental sensitivity to the weak dipole 
moments of the r are far from the MSSM predictions, at least three (two) 
orders of magnitude for the weak-magnetic (weak-electric) dipoles. 

Away from the Z peak both the electromagnetic and the weak dipole form 
factors are equally relevant but not enough to parameterize all the physical 
effects (in particular CP-violation). The case of f-pair production in high 
energy e + e _ colliders has been considered to illustrate this fact. Taking sev¬ 
eral CP-odd spin-observables based on the t and t polarization vectors, the 
different contributions from vertex and box corrections have been evaluated. 
There is no one loop contribution from the SM sector. It has been shown 
that, for the set of supersymmetric parameters that provides sizeable values 
of the t CP-violating dipoles, the SUSY box contributions happen to con¬ 
tribute with opposite sign and in a similar magnitude, yielding altogether a 
much smaller CP-violating observable effect. Another configuration has been 
shown for which the dipoles are smaller but the combined effect is larger. The 
same analysis has been performed using instead realistic observables based 
on the momenta of t and t decay products, with similar results. The chan¬ 
nels in which a bb , bd and id act as spin analyzers have been used under the 
assumption of standard CP-conserving decays and using the narrow width 
approximation. As expected, the leptons are the best t spin analyzers yielding 
the maximal values for the signal of CP-violation, r ~ 0.5 x 10 -3 . Nevertheless 
the statistics for such events is smaller. 
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A The one-loop tensor integrals 



Fig. A.l. Momentum convention for three-point-tensor integrals. 


The systematic way of dealing with one-loop integrals consists of reducing the 
tensor integrals to scalar ones. We employ the following set of tensor integrals 
(see Fig. A.l): 


C{0, At, nv}(PliP 2 , Mi, M 2 , M 3 ) 


167r 2 4 _ D r d D l {1, / M , 

i 11 J ( 2 n) D V{D 2 Vz 


(A.l) 


with 


T>i = (/ - pi) 2 - Ml + ie , 

X >2 — {l — P 2) 2 — Ml + ie , 

£>3 = l 2 - M 3 2 + ie . (A.2) 

and the orthogonal reduction (k± = p\ ±£>2) [76], 

C, = CtK, + Cfk., , (A.3) 

C^ w = C£ k + ^k +u + C 2 k_^k_ v + C'2 [k +ll k_ u + k +u k- fl ] + . (A.4) 


For us these definitions are very convenient, because the integral (A.l) is 
invariant under the combined replacements p\ p 2 and Mi <-»• M 2 , and we 
are dealing with equal mass final-state particles p\ = p\ = m 2 . Therefore 
Cf and C 2 + - are antisymmetric under these replacements while the other 
scalar integrals are symmetric, and the case Mi = M 2 automatically yields 

Cf=Ct= 0. 

For comparison we also list the tensor decomposition defined in [77]. 


L '{0, fi, 


(pi,p 2 ,M 3 ,M 1 ,M 2 ) 


16vr 2 4 _ D r d D l {1, / M , l^l u } 

i 11 J { 2 n) D V{D 2 Vz 


(A.5) 
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with 


Vi = l 2 — Ml + ie , 

T^2 — {l + Pi) 2 — -Wi + , 

©3 — {i + P 2) 2 — m'1 t ^ ■ 


and the tensor decomposition 


(A.6) 


c< D) =pi„c! D| +p 2 „cf > 


'/t — 1 1 (A-7) 

= QhuCqq ^ ^ + P2^P2 u C22 ^ + (Pl^iP2 u + ^ ■ 


(A.8) 


These scalar integrals are related to the ones obtained by orthogonal reduction 
in the following way: 


/'A- 0 ) _ /^»+ z^— 

L 1 — L 1 L 1 > 

cf ] = -c++cr, 

_ ^~i0 

L '00 — u 2) 

C{f ) =C 2 + + C 2 -+2C 2 + -, 

_/^»+ 1 o/^H— 

*~ / 22 — Ly 2 ~r l '2 ZL/ 2 I 

ri{D) _ (~t+ (~t— 

L-i 2 — o 2 — o 2 , 

with the arguments 

Cj D> = C, m (p 1 ,p 2 ,M 3 ,M I ,M 2 ), 


(A.9) 
(A-10) 
(A.11) 
(A.12) 
(A.13) 
(A.14) 




B Conventions for fields and couplings in the SM and the MSSM 

B.l The conventions 


We use the conventions of Ref. [78,79]. The covariant derivative acting on a 
SU(2 ) l weak doublet field with hypercharge Y is given by 


D 


v — 



(B.l) 


where r a are the usual Pauli matrices and the electric charge operator is 
Qf — /g + Y/2, with // = r 3 /2. The Z and photon fields are defined by 
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(B.2) 

(B.3) 


= W* cos 9 W - B M sin d w , 

Af, = sin 9 W + B M cos • 

The charged weak boson field is W± = 4^ (IT)] =|= iB'))). 

In the Standard Model (SM), there is only one Higgs doublet H with hy¬ 
percharge Y — 1. After spontaneous symmetry breaking (SSB), the physical 
Higgs field H° and the would-be-Goldstone bosons x and 0 ± are given by 


H = 


\M V + H ° + ix 'l 



(B.4) 


In the Minimal Supersymmetric Standard Model (MSSM) there are the same 
matter and gauge field multiplets as in the SM, each supplemented by su¬ 
perpartner fields to make up complete supersymmetry multiplets. The matter 
and gauge fields have the same quantum number assignments as in the SM. 
But in the MSSM there are two Higgs doublets with opposite hypercharges. 
Each forms a chiral supersymmetry nmltiplet and an SU(2) doublet. The La- 
grangian we use is defined in Ref. [80], especially in Eq. (1.45) and Eq. (1.54). 
We do not impose any reality constraint onto the parameters except for the 
reality of the Lagrangian. 

Spontaneous breakdown of the SU(2)xU(l) gauge symmetry leads to the ex¬ 
istence of five physical Higgs particles: two CP-even Higgs bosons h and H , a 
CP-odd or pseudoscalar Higgs boson A, and two charged Higgs particles H ± . 
They are grouped in two doublets, Hi = ir 2 &l and H 2 = 3> 2 , with opposite 
hypercharges (Y = =pl), where 


$i = 




(B.5) 


After SSB these are expressed in terms of the physical fields h , H, A, H ± and 
the would-be-Goldstone bosons G°, G ± by 


0+ = — H + sin (3 + G + cos (3 , (B.6) 

= —j= |ui + [(—hsina + Hcosa) + i(—Asin/3 + G°cos/3)]| , (B.7) 

4>2 = H + cos 0 + G + sin 0 , (B.8) 

02 = jv 2 + [(/icos a + H sin a) + i(Acos0 T G° sin 0)]} . (B.9) 
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Besides the four masses, two additional parameters are needed to describe 
the Higgs sector at tree-level: tan ft = v 2 /v\, the ratio of the two vacuum 
expectation values, and a mixing angle a in the CP-even sector. However, 
only two of these parameters are independent. Using M a and tan (3 as input 
parameters, the masses and the mixing angle a in the H, h sector read [79] 


M h,H ~ 2 + + e ) 


x 


IT 


\ 


1-4 


Mf Mf cos 2 2(3 + e(Mf sin 2 (3 + Mf cos 2 (3) 


(Mf + Mf + e) 2 


M h ± — Ma 


1 + 


Mi 


1 I/ 2 


W 




tan 2 a = tan 2/3 


Mf + Mf 


7r 


Mf — Mf + e/cos 2/3 ’ 2 <a< °’ 


which include the leading radiative correction 


(B.10) 

(B.ll) 

(B.12) 


e = ■ 

Here Gp 
squarks. 


3 G f 


ml 


v/2tt 2 sin 2 f3 


log 1 + 


m 2 / 


(B.13) 


is the Fermi constant and mq is the common mass scale for the 


The mass terms for the neutral gauginos and Higgsinos come out of the bilinear 
Higgs part of the superpotential (the /i term) and the soft-SUSY-breaking 
gaugino mass terms with masses M 2 and Mi for the SU(2) and U(l) gauginos 
A“ and X', respectively. The gaugino mass parameters are constrained by the 
GUT relation 


Mi = 5 -tM 2 ■ (B.14) 

3 ( ‘-\V 

Mixing terms arise from the minimal coupling terms between Higgs, Higgsino 
and gaugino fields, where the Higgs fields have been replaced by their vacuum 
expectation values. In terms of two-component spinors the mass terms add 
up to 


= --(-iV,-iA 3 ,^,^)U(-iA / ,—iA \r Hl ^k) T + h.c. , (B.15) 


with the symmetric mass matrix 
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0 


Y 


M 1 

0 


M 2 


\ 


—MzSw cos /3 MzCw cos /3 0 —/i 

^ M z sw sin (3 —M z cw sin (5 —fi Oy 


(B.16) 


We define the unitary diagonalization matrix N and the matrix N r , which is 
often more convenient, by the equations 


( 


Cw 


—sw 0 0 \ 


N*YN~ X 


N diag , N' = N ■ 


0 

0 


cw 0 0 
0 1 0 
0 0 l J 


(B.17) 


The matrix N is not unique. We could impose the requirement that N diag has 
only non-negative entries (which would force N to have an imaginary part if 
Y is not positive definite), but at least if the parameters in Y are real it is 
also possible to choose N real (which may lead to some negative entries of 
Ndiag)- However our formulas do not depend on this choice. In general the /j 
parameter may be complex and then the neutralino masses depend also on 
the phase of /j. The neutralino mass eigenstates are four Majorana spinors 
given by 


x] = {P L N jk YP R N* k )^ k , 



(B.18) 

(B.19) 


The mass terms for the charged gaugino and Higgsino have the same origin 
and, in terms of two-component spinors, read 


eg= 


/ 




v- 


x = 


0 X T 
X 0 

M 2 MwV% sin/3 
Mw \/2 cos ft n 




with the abbreviations 


(B.20) 
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(B.21) 


»/’/ = (- iX+ ,i 1 H,), 'Pj = (“iA ,<P Hl ) ■ 

We now define unitary diagonalization matrices U, V by the equation 


u* XV- 1 = M diag . (B.22) 

Again these matrices can be chosen to be real if X is real. Taking in general 
a complex /i parameter, the masses of the charginos are given by 

m~± 2 = \{Ml + \fi\ 2 + 2 Ml t [(Mf - \^\ 2 ) 2 + 4cos 2 2(3 (B.23) 

+4+ |^| 2 + 2M 2 Re(/i) sin 2/3)] ]1/2 }. (B.24) 

The chargino mass eigenstates are two Dirac spinors yb given by 


x+ = {p L v jk + p R u ; k )^ k , 



We abbreviate the charge conjugated fields as 


Xi 



(B.25) 

(B.26) 


(B.27) 


The mass terms for the scalar quarks come from the Yukawa couplings to the 
Higgs fields (which yield the corresponding quark masses) and the soft-SUSY- 
breaking squark mass terms and squark-Higgs interactions parameterized by 
rr iq L/R and A, respectively. Moreover there are mass and mixing terms from 
auxiliary field terms involving one or two Higgs bosons and two squark fields. 
The resulting mass matrix for the two squarks of a given flavor is (the de¬ 
pendence of the parameters rriq L/R and A on the generation/flavor has been 
dropped): 


M 2_f m L + m q rn* LR m q 
\ 'm L R'm q m\ + m\ 

where 

m 2 L = m 2 qL + cos 2(3 (l£ - Q f s 2 w )M 2 z , 
m \ = m \ R + cos 2/5 (QfSwWz » 



(B.28) 
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(B.29) 

(B.30) 



and, for up-type squarks: 


uilr = A — p* cot (3 , (B.31) 

while for down-type squarks: 

ttilr = A — n* tan /3 . (B.32) 

This hermitian mass matrix is diagonalized by a unitary matrix S, so we can 
write the squark mass eigenstates of flavor i as 


Qil — S u QiL + 5' 12 QiR , 

qi 2 = s\ 1 q iL + S\ 2 q iR . (B.33) 

In general vurr = \nriLR\e 1 ^ is a complex number and S may be written as 


e 19 ^ 2 cos 9i e ^ 2 sin d* 
y —e 1 ^/ 2 sin 9j e~^! 2 cos 9i 


(B.34) 


The mass terms of the sleptons arise in the same way as the squark masses. 
The main difference appears for the sneutrinos: there is only one sneutrino 
for every generation, iy, and hence there is no mixing. Moreover, we cannot 
add trilinear soft-breaking terms to shift the masses of the sneutrinos, whose 
value is given just by m 2 = rnf + 1/2 cos 2/3 M§. (See Eq. (B.29)). 

In general, the physical sfcrmion masses are given by 



m ) 



(m 2 R + m 2 L ) T 



m |) 2 + 4m 2 |m LR | 2 



(B.35) 


independent of the phase of uilr- 


B.2 Physical Phases in the MSSM with complex couplings 


As we do not constrain the parameters of the MSSM to be real, the following 
parameters may take complex values: the Yukawa couplings, /i, the soft pa¬ 
rameters m 2 2 , Mi, M 2 , M 3 and the A parameters. But not all combinations 
of phases in these parameters lead to different physical results because sev¬ 
eral phases can be absorbed by redefinitions of the fields. We now describe a 
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procedure to eliminate the unphysical phases. We assume the GUT relation 
between the M,, so they have one common phase. The only remaining phases 
will be chosen to be those of p, the A parameters and, for three generations, 
one phase for all the Yukawa couplings (the 5ckm)- 

Analogous to the Standard Model case, the Yukawa couplings can be changed 
by redefinitions of the quark superfields in such a way that there remains only 
one phase for three generations and only real couplings for less than three 
generations. 

Table B.l 

The charges rii for three U(l) symmetries that leave -Cmssm invariant. 


U(l) 

Mi 

A 

m\ 2 


Hi 

h 2 

QU 

QD 

LE 

6 

PQ 

0 

0 

-1 

-1 

1/2 

1/2 

-1/2 

-1/2 

-1/2 

0 

Ri 

-1 

-1 

0 

1 

0 

0 

1 

1 

1 

-1/2 


Following [12] we consider two U(l) transformations PQ and R\ that do not 
only transform the fields but also the parameters of the MSSM. In table B.l the 
charges n t of the various quantities are displayed with which the Lagrangian 
-Cmssm is invariant under the combined multiplications with e ian% . The first 
transformation is a Peccei-Quinn symmetry and R\ is an R symmetry that 
also transforms the 6 variables appearing in the arguments of the superfields 
in Cmssm- 

Since Cmssm is invariant under these transformations, so are the physical pre¬ 
dictions of the MSSM. If the parameters M*, A, m\ 2 and p have the phases 
(Pm, Ta, an d T\m we fi rs f a PPiy -Ri with the angle p>M, then PQ with the 
angle l fi m 2 12 and obtain for an arbitrary observable: 


0-(H, 1^1, \ M il l™l 2 l> 7W <PA, <PM,<Pml 2 ) 

= ^(M, \A\, \Mi\, \mj 2 \, ip ^ T Pm j Pa- Pm, 0, ^ m?2 ) (B.36) 

= v{\v\AA\,\Mi\,\m\ 2 \, (Pn + tpM-fPm^, <Pa-<Pm, 0, 0). 

So the physical predictions only depend on the absolute values of the param¬ 
eters and the phases 

4> A = arg(AM*), = arg- (B.37) 

One can replace these phases by another set where only the p and the A 
parameters are complex. In our choice the phase of A is traded for the phase 
of the off-diagonal term in the corresponding sfermion mixing matrix: rrr Lf{ = 
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Af — n* {cot,tan}/? (C.31,C.32). Relaxing universality for the soft-breaking 
terms, every Af contains a different CP-violating phase. 

B.3 Vertex factors 



P .AAAA* 7 = ieG(pi+p 2 )^ 

\ 

P2 X s 


N \ <|>2 



Fig. B.l. Generic Vertices. 

We employ the notation and conventions of the previous sections. The list of 
generic vertices is shown in Fig. B.l. 

B.3.1 Couplings in the SM 

Vertex 1: Coupling of electroweak gauge bosons and fermions: ie'y^V—A^f). 
lIf-V=-Q fx A = 0, (B.38) 
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(B.39) 


Zff : V 


W ± ff-.V 


v f A _ a f 
2 S\yCw 2 SyyCw 


2\[2sw 


(B.40) 


with v f = (it - 2slrQf), a f = /•{. 

Vertex 2: Coupling of gluons and quarks: —ieg s 'y p T a . 

Vertex 3: Coupling of two fermions and one Higgs boson (7 Iff): ie(S— -P 75 ). 


H°ff: S = ~n f , P = 0, (B.41) 

xff ■ S = 0, P = -2ilifi f , (B.42) 

0 ± /7 : -5 = V2l/ ifif - / 17 ] , P = -V 21 I \g f + /^] , (B.43) 

rn f rrif 

with iif = - - = - - •' r . 

2 S\yAI\y 2swC\yMz 

Vertex 4: Coupling of one Higgs boson and two gauge bosons: i eKg pu . 


H°ZZ : K = —- , (B.44) 

SwCw 

H°W ± W T :K = — , (B.45) 

% 

(j) ± W T 7 : K = M w , (B.46) 

cj) ± W T Z : K = —M z s w ■ (B.47) 

The rest of the couplings are K = 0. 

Vertex 5: Coupling of one gauge boson and two Higgs bosons (V0i0 2 ) : 


ieG(pi + p 2 ) P - 



7<TV _a 

:G = -A , 

(B.48) 

Z(j) x (j)- x 

x cos 2 9 W 

2 SwC\y 

(B.49) 

Z X H° 

■ G = 2 S L’ 

(B.50) 

O 

1 

-e- 

■ G “ 2 Sw ’ 

(B.51) 

w x r x x 

: G = 

2 sw 

(B.52) 


Interchanging the two Higgs bosons causes the coupling constant to change 
sign. 

Vertex 6: Coupling of three gauge bosons (outgoing momenta): 

i eJ{g^(k x - k°) p + g up (k~ x - k x ) p + g pp (k° - k~ x ) v ) . 

X W x W~ x : J = -A , (B.53) 
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(B.54) 


ZW X W~ X :J = -. 


B.3.2 Couplings in the MSSM 


Vertex 1: Coupling of two neutralinos or charginos and one gauge boson: 

ie^V - A 75 ). 

The fermion flow direction in our Feynman graphs is fixed by the outgoing 
fermion antifermion pair. The relation between the gauge boson vertex fac¬ 
tors for the two different fermion flow directions is obtained by substituting 
spinors by charge conjugated spinors in the interaction term: 


V,x + Y(Pl9l + Pr9r)x + = VnX 7 m (~Pl9r ~ Pr9l)x • (B.55) 


Taking into account also the symmetry factor S = 2 for the neutralino 
coupling and S — 1 for the charginos, the boson vertex factors V = (gi, + 
9r)/2 and A = (g L - g R )/ 2 are given by: 


ZXiXl : 9l = 


-(K 4 K 


/* 
j 4 


2 S\yCw 

9r = —-— {N'CN' 
yK 2s w c w J 


K 3 n;i) , 


Zxtxi ■9 l = 

9R = 
z Xj Xk '■ 9l — 
9r = 


1 


SwCw 

1 


Sw c W 

1 r 


SwCw 

1 


SwCw 



- *w) v; 2 v )2 + 4,i4*Vii] , 

— j Un'J}, — 'uW: 1 

Uk 2 Uj2 + 4 u 2 1 r"* , 

v ki v n + 4UW).] , 


and 


(B.56) 

(B.57) 

(B.58) 

(B.59) 

(B.60) 

(B.61) 


7 XkXj ■ 9 l = 9r = Sjk , (B.62) 

7 Xj Xk '■ 9 l = 9r = h jk • (B.63) 


Vertex 2: There is no genuine supersymmetric vertex of this kind. 

Vertex 3: There are three couplings of the kind ie(S' — P 75 ) : 

(Some of the couplings are more easily written in terms of 9l,r = S ±P). 
• Coupling of one neutralino or chargino to fermions and scalar fermions. 
The couplings of neutralinos and charginos to quarks and scalar quarks 
are given by 


P-xqq = e Q ] x{PL9L + Pr9r)9 + eq(P R g* L + P L 9 R )XQ • (B.64) 


down-type quarks (i = 1,3,5): 
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Qi,kXjQi ■ 9l 


= -V2 


QiN'jl + (I l 3 - Q z s 2 w ) - n'* 2 si, 

SwCw ' 


+ 


m, 


<n 


-N^Si 


'2M w s w cos (5 i3 k2 
9r = — V2 


Qi N jl + (~QiSw) - N f2)Sk 2 

SwCw 


_ mqi ]\T' Qi 

2M w s w cos(3 j3 kl 


9(i+i),kXj 9i ■ 9l 


Sw 


m 

T r* C*+l 1 Qi+l ,/■* c*+l 

V jl°kl ^ /K nr . a V j2 D k2 

\/2M w snip 


9r — — 
Sw 


m, 


Qi 


u i 2 si +1 


_V2M w cos(3 j2 kl 
up-type quarks (i = 2,4,6): 

= -V2 


qlkXjQi ■ 9 l 


QiNjl + (II - QiS 2 w )- N» SI 

SwCw ' 


kl 


m, 


ii 


2 M w sw sin /3 


N'* < 7 * 

Jv j'4°fc2 


9R 


= -V2 


— (QiNj 1 + (— QiSw) - N j2 ) S l k2 

\ SwCw / 


m, 


<u 


2M w sw sin f3 


vl-1),kX-Ui'■ 9L 


s w 


N'jA 


777 

TT* Qi— 1 | Qi—1 jj * Qi— 1 

+ V2M w c OS l3 U * bt2 


9r — — 
% 


m, 


n 


y/2Mw sin f3 


vpsi i 1 


(B.65) 


(B.66) 

(B.67) 

(B.68) 


(B.69) 


(B.70) 

(B.71) 

(B.72) 


The couplings of neutralinos and charginos to leptons and scalar leptons 
are given analogously, performing the following substitutions: 


i = 1,3,5 :qn - 

h j 

ni 

^jk 

-j. Q l 

Qi = - 1 , 

qi2 ~ 

+ h , 




i = 2,4,6 :q a - 

-»■ v\ , 

Qi 
^jk 

djk , 

Q* = 0 , . 


, = m , 


"Qi 


qi 2 does not exist . 


Coupling of one gluino to a quark and a scalar quark. 

The interaction between gluinos, quarks and squarks is described by the 
terms 


£ = q ] eg a (P L g L + P R g R )yg + eq(P R g* L + P L g %)y , 


(B.73) 
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yielding the vertex factors 


ie( p L g L + Pr9r)y^ ie ( p R9*L + p l9r) y- (B.74) 

In our calculations the Gell-Mann matrices appear only in the combination 
8 /\ a \ a \ A 

Eh- = C 2 (F)S ab = -Sab. (B.75) 

a=l \ 4 / AB 6 

The couplings are 

■ egL = -V2g s Sn , (B.76) 

eg R = +V2g s S i2 . (B.77) 

• Coupling of two fermions and one Higgs boson (Tiff). 


Huu 

S = — ji u sin a / sin f3 , 

P = o, 

(B.78) 

Hdd 

S = —/irfcoscc/cos/3 , 

P = 0, 

(B.79) 

huu 

S = — p, u cos a / sin (3 , 

P = 0, 

(B.80) 

hdd 

S = gd sin a / cos /3 , 

P = 0, 

(B.81) 

Auu 

5 = 0, 

P = —\g u cot (3 , 

(B.82) 

Add 

5 = 0, 

P = —i/irftan/3 , 

(B.83) 

H + ud 

5 = —p(/x n cot /3 + /irftan/3) 
V2 

1 

P = —j= (n u cot (3 — Hd tan (3) 
V2 

? 

5 

(B.84) 


G°ff '■ 5 = 0, P = -‘hllji, . (B.85) 

G ± f'S- S=s/2I&P — + Hf>] . (B.86) 


with jif = rrif /2s\yMw = mf/2swCwMz- For the vertices corresponding to 
the hermitian conjugated fields we have to replace (. S,P ) by (S *, —P*). 
Vertex 4: Coupling of one Higgs boson and two gauge bosons: i eKg^ u . 

The couplings of a neutral Higgs to two Z bosons are: 

My 

hZZ : K = -— sin(/3 - a) , (B.87) 

SwCw 

My 

HZZ : K = -— cos(/5 - a) . (B.88) 

SwCw 

And the only nonzero couplings of a charged Higgs to neutral gauge and W 
bosons are: 


G^H^y : K = M w , 
G^^Z : K = —M z s w . 
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(B.89) 

(B.90) 




The rest of the couplings are K = 0. 

Vertex 5: There are two types of couplings Zcp i0 2 of the kind ieG(pi + p 2 )// 


Coupling of two Higgs bosons and one neutral gauge boson. 

i sin(/3 — a) 


ZAH:G= 
ZAh:G= 
ZG°H : G = 
ZG°h : G = 
ZH X H~ X :G = -X 
ZG X G~ X :G = -A 


2S]yCw 
i cos (/3 — a) 
2 SyyCw 
i cos (/3 — a) 
2s\yC\y 
i sin(/3 — a) 

2s\yCw 
cos 2 6 W 


2 S]yCw 
cos 28 w 
2S\yC\y 
^H X H~ X : G = -A , 
7 G a G- a : G = —X . 


(B.91) 

(B.92) 

(B.93) 

(B.94) 

(B.95) 

(B.96) 

(B.97) 

(B.98) 


Interchanging the (pi and 0 2 causes the coupling constant to change sign. 
Coupling of two scalar fermions and one neutral gauge boson. 


'yQi,jQi,k • ^ Qi^jk i 

ZQijQi,k '■ G = 


(Is - QiSw)S7i 1 2 3 * S l i - QistrSZSl 2 


SwCw 

Vertex 6: There is no genuine supersymmetric vertex of this kind. 


(B.99) 

(B.100) 
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